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SUMMARY 

The  regional  and  teleseismic  signals  from  a  set  of  short  period  stations  of  the  Israel 
Seismic  Network  were  investigated  with  the  primary  objective  of  seismic  source 
discrimination.  The  regional  database  includes  212  ripple-fired  and  single  shot  quarry  blasts, 
underwater  explosions  and  earthquakes  of  local  magnitudes  Ml=1. 0-3.1,  at  distances  5-250  km 
from  various  areas  in  Israel,  Jordan  and  the  Mediterranean  Sea.  Teleseismic  data  consist  of 
68  earthquakes  and  nuclear  tests  from  Eurasia  (mb=5.2-r6.5).  Also  included  in  the  study  are 
NORESS  recordings,  comprising  36  events,  almost  completely  intersecting  with  ISN  dataset. 

The  multichannel  seismograms  and  their  smoothed  spectra,  calculated  for  the  whole 
signal,  demonstrate  several  specific  features  (S-wave  group  velocity,  coherency  of  spectral 
shapes,  spectral  energy  distribution)  essentially  different  for  regional  earthquakes,  quarry  blasts 
and  underwater  explosions.  These  features  were  used  to  construct  several  new  efficient  scalar 
and  vector  statistics  calculated  in  specific  frequency  bands:  spectral  semblance  (SS),  energy 
spectral  ratio  (SR)  of  low  and  high  frequencies  and  spectrograms,  as  well  as  velogram 
parameter  C  (VPC)  estimated  in  the  time  domain. 

Owing  to  dominant  low-frequency,  low-velocity  surface  waves  and  source  spectral 
modulation  in  seismograms,  quarry  blasts  and  underwater  explosions  showed,  as  a  rule,  higher 
SR  and  SS  values,  and  lower  VPC  as  compared  to  local  earthquakes.  Different  frequency  bands 
were  tested:  the  ratio  (1-3  Hz)/(6-8  Hz)  and  semblance  (1-12  Hz)  provided  the  best 
discriminant  performance,  with  a  success  rate  of  97%  for  the  whole  database. 

The  majority  voting  as  well  as  automatic  classifiers:  the  Fisher  Linear  Discrimination 
Function  (LDF)  and  the  non-linear  Artificial  Neural  Network  (ANN)  applied  to  combinations 
of  the  above  mentioned  discriminants,  were  tested  on  quarry  blasts  and  earthquakes  in  three 
separate  regions  of  Israel  and  Jordan.  The  performance  of  the  different  classification  schemes, 
determined  by  the  {classifier,  parameter  vector}  pair,  was  compared  on  the  basis  of  error 
probability  estimation  using  the  leave-one-out  cross-validation  procedure.  The  best  result,  with 
an  error  probability  of  0.006,  was  achieved  when  ANN  was  applied  to  the  5-band  spectrogram 
in  the  (1-1 1  Hz)  frequency  range,  averaged  over  an  ISN  subnetwork. 

The  SS,  the  SR  and  the  ratio  of  P/P-coda  spectral  maxima  were  effective  for 
classification  of  teleseismic  events.  Contrary  to  observations  of  local  data,  SS  and  SR  values 
calculated  from  ISN  seismograms  in  time  windows  containing  P  and  P-coda  waves,  are  higher 
for  earthquakes.  The  best  performance  is  provided  in  the  frequency  ranges  (0.6-1  Hz)/(l-3  Hz) 
for  SR  and  (0.6-2  Hz)  for  the  SS.  Joint  application  of  the  two  discriminants  showed  95% 
separation  between  the  two  populations.  None  of  the  individual  statistics  resolved  nuclear  tests 
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and  earthquakes  completely,  but  majority  voting  on  some  combinations  provided  total 
separation.  The  leave-one-out  cross-validation  test  revealed  the  LDF  as  the  best  automatic 
classifier  (only  one  misidentification  for  the  dataset  without  deep  earthquakes). 

Teleseismic  ISN  and  NORESS  array  discrimination  performances  were  compared.  The 
NORESS  records  showed  the  same  (as  the  ISN)  relation  between  SR  values  for  earthquakes 
and  explosions,  whereas  the  character  of  SS  was  reversed.  By  applying  (0. 6-0.8  Hz)/(0.6-10 
Hz)  ratio,  we  obtained  full  separation  of  the  two  populations.  In  general,  SR  perfoimance 
appears  to  be  better  for  NORESS  recordings,  whereas  SS  is  better  for  the  ISN  data. 

We  tested  array  methods,  applied  to  the  different  ISN  configurations  and  investigated 
teleseismic  detectability  and  slowness  estimation.  The  standard  beam  (SBT)  and  adaptive 
beamforming  technique  (ABT)  were  compared  for  ISN  subnetworks  containing  7-16  stations. 
Theoretical  analysis  and  experimental  results  yielded  that  only  ABT  is  capable  of  improving 
the  detectability,  providing  the  best  SNR  at  the  output.  Estimation  of  slowness  using  the 
adaptive  optimal  detector,  maximum-likelihood  onset  estimator  and  phase  correlation  analysis, 
provided  deviation  of  slowness  as  <0.003  sec/km  and  azimuth  <0.5°  for  two  small  events 
(mb<5.4). 

Some  physical  phenomena  and  experimental  properties  related  to  local  explosive  sources 
were  analyzed  based  on  ISN  recordings.  Unique  spectral  features  of  seismic  waves  from  the 
deep  experimental  detonations  in  the  Dead  Sea,  caused  by  source  (bubble  oscillation)  and  path 
propagation  (reverberations)  effects,  were  observed  and  interpreted.  The  two  effects  produce  a 
distinct  spectral  modulation  by  two  series  of  harmonics  at  frequencies  corresponding  well  to 
those  estimated  from  the  source  parameters  and  the  physical  properties  of  the  Dead  Sea  water. 

Seismic  efficiency  of  local  explosions  and  peak  amplitude  attenuation  were 
experimentally  investigated.  A  clear  correlation  between  charge  weight  and  magnitude  is 
observed  for  instantaneous  quarry  shots  and  single  underwater  explosions  (unlike  ripple-fired 
blasts  with  specific  delays  and  geometry).  A  dependence  of  peak  particle  velocity  on  the  scaled 
distance  was  evaluated.  A  high  coupling  for  Dead  Sea  UWE  is  also  observed:  an  0.3  ton  charge 
corresponds  to  Ml=3.1,  exceeding  magnitude  values  of  quarry  blasts  of  several  dozens  tons. 
This  observation  should  be  considered  when  planning  calibration  explosions,  currently  lacking 
in  the  region,  and  necessary  for  the  improvement  of  CTBT  monitoring.  According  to  the 
equations  developed,  a  Dead  Sea  UWE  of  5,000  kg  will  “produce”  Ml  >  4.0  and  might  be 
observed  at  distances  up  to  2600  km,  i.e.  might  be  recorded  by  most  of  seismic  stations  and 
local  networks  operating  in  the  Eastern  Mediterranean  and  Middle  East. 
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1.  INTRODUCTION 


Seismological  means  of  detecting  and  identifying  nuclear  tests  are  the  most  powerful 
methods,  but  still  need  further  improvement.  Under  the  Comprehensive  Test  Band  Treaty 
(CTBT),  monitoring  of  seismic  events  as  small  as  mi,=2.5  may  become  a  crucial  issue  (Wallace 
et  al.,  1992),  especially  when  considering  evasion  scenarios  such  as  decoupling  (Murphy  et  al., 
1987).  The  weak  events,  including  those  resulting  from  human  activities  such  as  construction 
and  mining,  must  be  discriminated  from  nuclear  explosions  (Stump  et  al.,  1994). 

The  wavefields  of  these  small  events  are  only  well  recorded  at  regional  distances,  thus 
attracting  interest  to  regional  seismic  networks.  The  ISN  is  one  of  the  best  in  the  Middle  East, 
where  international  studies  concerning  the  CTBT  have  just  begun.  Physical  factors  affecting 
regional  discrimination  of  man-made  and  natural  seismic  events  (generation  of  dominant  low- 
frequency  surface  waves  by  quarry  blasts,  associated  with  subsurface  sediments;  differences  in 
spectral  content,  kinematic  parameters  of  seismograms  and  spatial  coherency  of  spectral 
shapes)  were  investigated  by  Kafka  (1990),  Bath  (1975),  Alexander  etal.  (1995),  Gitterman 
and  van  Eck  (1993)  and  others.  The  physical  phenomena  observed  at  a  regional  dense  seismic 
network,  provide  a  basis  for  development  of  efficient  classification  algorithms,  optimally 
accounting  for  the  multichannel  character  of  the  system. 

A  successful  Integrative  Approach  using  majority  voting,  artificial  neural  network,  linear 
Fisher  and  other  multivariate  classifiers  applied  to  combinations  of  discriminants,  was 
demonstrated  (e.g.  Wuster,  1993;  Fully,  1995  and  Dowla,  1995)  to  further  enhance  the 
resolving  power  of  discrimination.  The  optimal  parameter  vectors  and  best  classifier  algorithms 
can  be  elaborated  on  the  basis  of  error  probability  estimation  (Tsvang  et  al.,  1993)  through  a 
leave-one-out  cross-validation  procedure  (Fukunaga  and  Hummels,  1989). 

The  ISN  potential  to  contribute  to  teleseismic  monitoring  has  not  been  evaluated. 
Enhancement  of  the  ISN  detectability  can  be  achieved  through  the  application  of  the  array 
technique  to  its  subnetworks.  New  adaptive  optimal  detector  and  beamforming  (Kushnir  et  al., 
1990)  together  with  phase-correlation  technique  (Gelchinsky  et  al.,  1985)  were  applied  to  the 
data  yielding  additional  noise  suppression  and  accurate  slowness  estimates  for  precise  location 
of  weak  teleseismic  events. 

A  number  of  CTBT  studies  claim  lack  of  calibration  explosions  in  the  Middle  East  region 
and,  in  particular,  for  areas  along  the  Dead  Sea  transform  (Baumgardt,  1996).  Design  and 
execution  of  such  experiments  require  empirical  assessment  of  explosion  effectivity  and 
attenuation  of  seismic  waves  for  regional  geological  and  geotechnical  settings. 


2.  DATABASE. 


2.1  Regional  Events. 

Regional  data  include  seismograms  from  212  quarry  blasts,  underwater  explosions 
(UWE)  and  local  earthquakes  of  magnitudes  Ml=1. 0-2.8  and  depths  H=0-28km  from 
geologically  different  areas  in  Israel,  Jordan  and  the  Mediterranean  Sea.  The  events  were 
recorded  at  distances  of  5-250  km  by  the  vertical  short-period  stations  of  the  Israel  Seismic 
Network  (50  Hz  sampling  rate,  0.2-12.5  Hz  bandpass  filtering).  We  examined  all  the  traces  and 
eliminated  records  that  appeared  to  be  affected  by  glitches,  spikes,  noise  bursts,  or 
malfunctions.  Most  explosions  were  supplied  by  ground-truth  information  acquired  from 
quarry  operators.  Basic  parameters  of  the  selected  events  and  information  about  the  quarries  are 
presented  in  Tables  1  and  2  and  the  locations  are  shown  on  a  map  (Fig.  1).  Complete 
information,  regarding  earthquake  parameters,  blast  groimd-truth  information  and  event 
numbers  used  in  this  report,  is  available  from  the  Scientific  Report  No.  1  (Gitterman  et  al., 
1996b);  data  about  Jordanian  quarries  are  taken  from  the  USGS  report  (Wilburn  et  al.,  1996). 


Table  1.  Local  database  used  in  the  study 


# 

Region,  dataset 

Number  of 
Earthquakes 

Number  of 
Explosions 

Comments 

1 

Galilee,  Northern  Israel 

30 

39 

Earthquake  depth  H  =  0-21  km 
Quarry  blasts:  W=2000-14000  kg 

2 

Southern  Dead  Sea  and 
Negev  desert,  Israel 

16 

26 

H  =  0  -  24  km 

Quany  blasts  W= 1700- 105 00  kg 

Gilad  region,  Jordan, 
Israel-Jordan  border 

19 

15 

H  =  2  - 18  km 

Quarry  blasts  -  no  GTI 

1 

Tyre  region, 
Mediterranean  sea 

8 

16 

H  =  6  - 17  km 

UWE  -  unauthorized  fishing,  no 
Ground  Truth  Information 

5 

Dead  Sea  basin,  seismic 
refraction  profile 

- 

28 

UWE  for  seismic  prospecting 
depth  h  =  70  m;  W=1 6-304  kg 

6 

Jordan/Saudi  Arabian 
border 

- 

5 

Jordanian  quarry  blasts 

W=  15000-17000  kg 

7 

Quarries  Revaya, 

Mt.  Nitzim 

- 

10 

Single  quarry  blasts  (no  delays) 

W  =  900  -  2800  kg 
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Fig.  1 .  Locations  of  ISN  stations,  quarries  and  local  seismic  events  used  in  the  study. 


Table  2.  Information  about  open-pit  quarries  presented  in  the  study. 


II 

r  '  NAME 

CODE 

LAT. 

LONG. 

PRODUCT 

ISF 

LAEL 

■ 

Aroar 

31.125 

35.003 

phosphates 

2 

Carmel 

nshr 

32.758 

35.040 

dolomites  (construction,  roads) 

H 

Dragot 

drgt 

31.309 

35.078 

phosphates 

1 

Golani 

goln 

32.782 

35.397 

aggregate,  blocks,  beton 

5 

Hanaton 

hntn 

32.773 

35.241 

dolomites  (construction,  roads) 

6 

Kadarim 

kdar 

32.913 

35.465 

limestone,  dolomite  (construction, 
roads) 

7 

Mt.Nitzim 

ntzm 

31.184 

35.051 

calcium  carbonates 

8 

Oceana 

31.208 

35.385 

salt  aggregates 

9 

PhosphArad 

31.131 

35.210 

phosphates 

PhosphZin 

^.851 

35.082 

phosphates 

Poria 

pori 

32.764 

35.553 

basalt,  aggregate,  asphalt,  beton 

H 

Revaya 

rvya 

32.439 

35.480 

limestone  (construction,  roads) 

H 

Shefaram 

shfr 

32.797 

aggregate,  asphalt,  beton 

H 

Tamra 

tmra 

32.858 

35.220 

carbonate,  chalk,  cement 

Vered 

verd 

32.479 

35.069 

dolomites 

JOR 

DAN 

II 

Jarash 

32.28 

rail 

2 

Wadi  Mujib 

31.47 

35.57 

gypsum,  sulphur 

Gilad 

32.46 

35.86 

stone 

Some  of  the  collected  data  are  included  in  the  Ground-Truth  Database  (GTDB), 
developed  by  Multimax  and  presented  on  the  Web  site  (see  Fig.  2): 

http://esl.multimax.com/~gtdb/galilee/index.html 
(Grant  et  al.,  1997).  The  data  were  used  by  other  researchers  in  a  discrimination  study  (Walter 
etal.,  1997)  and  in  an  estimation  of  chemical  explosions  effectivity  (Khalturin  and  Richards, 
1997) 
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«  GTDB:  Gn/GalUee  Dataset  Home  Page 


select  earthguakes  or  quarry  blasts  from  the  maps  below 


r  Earthquakes  I  Quarry  Blasts  1 


How  to  Get  the  Galilee  Data  Set 
About  the  Israeli  Seismic  Network  ^  About  the  Galilee  Dataset 

Galilee  Home  I  News  I  Credits  I  References 


lori@otter.mmiwest.com 


Fig.  2.  The  Web  page  prepared  by  Multimax  with  Israel  local  events  selected  and 
checked  in  the  project. 
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Spectral  computations,  necessary  for  estimating  spectral  discriminants  used  in  the  study, 
are  made  for  a  time  window  of  about  20-30  sec,  which  is  the  same  for  all  stations  for  a  fixed 
event,  and  includes  the  whole  signal  up  to  the  end  of  the  coda  waves.  This  window  was 
selected  to  avoid  picking  out  separate  wave  phases  at  short  distances  and  accumulated 
information  regarding  source  features  being  kept  in  all  wave  forms,  thus  enhancing  the 
resolving  power  of  the  discriminants  considered.  The  spectra  were  instrument-corrected  and 

smoothed  m  a  fixed  0.5  Hz  moving  window  to  provide  equivalent  spectral  resolution  for 
different  stations. 

2.2  Teleseismic  Events 

The  teleseismic  database  consists  of  29  nuclear  tests  from  the  test  sites  in  China  (Lop 
Nor),  Kazakhstan  (Semipalatinsk),  and  Russia  (W.  Siberia,  N.Zemlya),  12  deep  (H>50km) 
earthquakes  from  Afghanistan  and  27  shallow  earthquakes  from  China  and  Tajikistan.  These 
events  were  recorded  by  the  ISN  and  their  magnitudes  are  in  the  range  mb=5.2-^6.5.  (For 
several  Chinese  nuclear  explosions,  we  used  in  addition  recordings  of  the  Israel  broad-band 
station  BGIO).  In  parallel,  during  the  visit  to  NORSAR  we  collected  and  processed  NORESS 
seismograms  of  20  explosions  and  14  earthquakes  and  two  additional  earthquakes  (non- 
readable  on  ISN  tapes).  Other  events  available  to  the  ISN  were  inaccessible  in  the  NORESS 
database.  Detailed  source  information  is  presented  in  Tables  3  and  4;  event  locations  are  shown 
on  Fig.  3.  Location,  origin  time,  magnitude  (m^)  and  focal  depth  are  taken  from  PDE  (NEIC) 
bulletins,  distance  and  azimuth  are  calculated  relative  to  the  station  BGIO  of  the  ISN. 

The  records  (2-2.5  minutes  in  length)  contain  only  P-wave  and  P-coda  waves  owing  to 
the  restncted  recording  time  of  the  acquisition  system.  A  time  window  about  13-15  sec 
comprising  the  P-signal,  was  used  for  the  spectral  analysis,  the  spectra  were  smoothed  in 
0.2  Hz  window.  ISN  spectra  were  instrument-corrected. 

3.  ANALYSIS  OF  LOCAL  EXPLOSIVE  SOURCES. 

3.1  Some  Spectral  F eatures  of  Underwater  Explosions  in  The  Dead  Sea. 

A  series  of  single-shot  underwater  explosions  with  charges  of  W= 16-3 04  kg  was 
detonated  in  the  Dead  Sea  at  a  depth  of  d=70  m.  For  the  analysis  we  selected  a  representative 
dataset  of  28  explosions  of  four  different  charge  sizes,  distributed  uniformly  along  a  line  in  the 
middle  of  the  northern  basin  of  the  Dead  Sea  (see  map  in  Fig.  1).  Ground-truth  information, 
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Test  sites  &  earthquakes 


Fig.  3.  Locations  of  nuclear  test  sites  and  teleseismic  events  used  in  the  study.  Deep 

(Afghan)  earthquakes  are  indicated  by  darker  circles.  Israel  Seismic  Network  is 
indicated  by  the  station  BGIO  (Bar-Giyyora)  placed  in  the  middle  of  ISN. 
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including  accurate  coordinates,  charge  weight,  shot  depth  and  local  time  of  explosions  is 
presented  in  Gitterman  et  al.  (1996). 

Detonations  in  deep  water  produce  seismic  waves  bearing  spectral  features  caused  by  the 
source  (bubble  oscillation)  and  path  propagation  effect  (reverberations).  The  bubbling  effect 
produces  a  low-frequency  azimuth-invariant  spectral  modulation  as  the  complete  harmonic 
series  of  spectral  maxima  with  a  fundamental  frequency  fjj,,  which  is  characterized  by  the 
detonation  depth,  d  (m),  and  the  explosive  yield,  W  (kg  of  TNT)  (see  e.g.  Willis,  1963): 

fi,  =  (d+10)“/(2.iw''’);f„,  =  n.4,  n=l,2,...  (3.1) 

The  bubbling  actually  affects  all  types  of  recorded  seismic  waves  throughout  the  entire 
wavetrain  shown  in  Fig.  4,  but  is  better  manifested  in  spectra  of  the  initial  body-wave  portion 
{P  and  P-coda).  A  complete  harmonic  series  (3,  6,  and  9  Hz)  resulting  from  the  bubbling  effect 
is  presented  in  Fig.  5a  for  the  largest  explosion  EXl  with  charge  weight  W=304  kg.  Spectra  of 
the  whole  signal  are  much  more  complicated  and  mixed  with  reverberations  maxima  (Fig.  5b), 
however  low-frequency  spectral  modulation  minima  (4.5  and  7.5  Hz),  corresponding  to  the 
bubbling  interference,  are  sharpened  and  expressed  more  distinctly. 

Using  ground  truth  information  and  taking  into  account  the  high  density  of  salt  water  in 
the  Dead  Sea  i.e.  p«l.ll  gr/cm^  compared  to  the  usual  p»1.03  gr/cm^  (Prof.  Z.  Ben-Avraham, 
oral  communication)  we  estimated  from  Eq.  3.1  the  fundamental  bubbling  frequency  for 
explosion  EXl  as  fb«2.9  Hz.  This  corresponds  well  to  the  observed  value  ~3.0  Hz  presented  in 
Fig.  5a.  (The  usual  density  yields  an  estimate  4  w2.7  Hz). 

In  contrast  to  the  complete  harmonic  series  caused  by  the  bubbling  effect,  water 
reverberations  produce  two  odd  harmonic  series  with  fundamentals  fn  and  4  (Weinstein,  1968): 

fH=  Vo/4H,  4b  =  (2n-l)*  4,  n  =  1, 2,...  (3.2) 

fd=Vo/4d,  4b  =  (2n-l)*  4,  n=l,2,...  (3.3) 

where  H  is  the  water  depth  in  meters  and  d  is  the  detonation  depth  in  meters. 

For  the  P  velocity  in  the  Dead  Sea  Vo~l 770  m/sec,  and  for  an  average  depth  of 
H  «  310  m  (Neev  and  Hall,  1979),  we  estimated  from  Eqs.  2  and  3,  fH»1.4  Hz,  and  4«6.3  Hz. 
The  spectral  modulation  caused  by  reverberations  is  better  manifested  in  the  spectra  of  the 
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Fig.  4.  Seismopams  of  Dead  Sea  underwater  explosions:  (a)  the  largest  explosion 
EXl  with  charp  weight  W=304  kg;  arrows  show  a  coda  window  used  in  the 
analysis  (see  Fig.  5c);  (b)  a  small  explosion  EX13  showing  some  earthquake 
wave  features  (e.g.  the  largest  amplitudes  of  the  S/Rayleigh  arrivals  at  each 
station).  Amplitudes  at  this  and  following  figures  are  not  scaled. 
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coda  window  shown  on  Fig.  4a.  The  spectra  on  Fig.  5c, d  demonstrate  distinct  maximum 
frequencies,  fitting  the  estimated  fn  ,  odd  harmonics  Sfn  and  Tfn  and  a  less  visible  harmonic 
Sfn-  A  noticeable  maximum  at  about  6  Hz  contains  supposedly  joint  contributions  of  the  close 
Sfn  (7  Hz),  fjj  (6.3  Hz)  and  24  (5.8  Hz)  harmonics. 

Correct  interpretation  of  spectral  features  of  recorded  seismic  waves  is  useful  in 
evaluating  energy  distribution  in  the  frequency  domain  depending  on  different  effects,  source 
and  medium  parameters  and,  thus,  in  estimating  spectral  discriminants  performance. 

Simple  comparison  of  the  seismograms  of  explosions  (Fig.  4a)  and  earthquakes  (Fig.  6a) 
shows  differences  that  are  due  to  features  of  the  source.  For  the  “typical”  earthquake 
seismograms  the  ■S'/Rayleigh  waves  can  usually  be  identified  and  have  the  largest  amplitudes  at 
each  station.  For  “typical”  explosions  (especially  for  UWE),  the  largest  amplitudes  are  usually 
in  the  P-group  with  57Rayleigh  onsets  being  difficult  to  identify.  However,  there  are  some 
events  which  demonstrate  anomalous  amplitude  characteristics;  for  example,  a  small  explosion 
EX13  shows  some  earthquake  wave  features,  i.e.  the  largest  amplitudes  of  the  ^/Rayleigh 
arrivals  at  each  station  (Fig.  4b)  and  explosion-like  wave  features  for  an  earthquake  QS12  can 
be  observed  (Fig.  6b).  For  reliable  identification  of  seismic  events  in  the  region,  therefore, 
simple  visual  analysis  of  seismograms  is  not  sufficient  and  require  additional  spectral  analysis 

3.2  Experimental  Relationships  of  Explosion  Effectivity  and  Attenuation. 

There  are  well  established  correlations  between  the  magnitude  and  yield  of 
underground  nuclear  explosions  for  different  test  sites  (e.g.  Ringdal  et  al.,  1992): 

mb  =  4.45+0.75*log,oY(kton)  (3.4) 

We  have  checked  similar  relations  for  different  types  of  blast  sources  presented  in  our 
database.  Most  studies  of  ripple-fired  quarry  blasts  show  no  correlation  between  total  charge 
weight  of  a  blast  and  magnitude  (or  PGV  -  peak  ground  velocity),  one  example  is  the 
investigation  by  Jarpe  etal.  (1996),  who  studied  production  and  waste  shots  of  10-180  tons, 
with  delays  of  mostly  65  and  100  msec  between  30-100  short  rows  of  charges  oriented 
diagonally  to  the  bench,  at  the  Carlin  open-pit  gold  mine.  However,  in  the  case  of  big  Kursk 
mining  blasts  of  200-1000  tons  with  delays  of  mainly  20  and  35  msec  between  three  long  rows 
of  charges  oriented  along  the  bench,  a  clear  dependence  of  magnitude  on  total  charge  was 
observed  (Leith  etal.,  1997).  Following  a  simple  analysis  of  superposition  and  interference  of 
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Fig.  6.  Examples  of  Dead  Sea  earthquakes:  (a)  “typical”  seismograms  of  QS3; 
(b)  explosion-like  wave  features  for  QS12  (the  largest  amplitudes  in  the  P-group  for 
the  five  closest  stations). 
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seismic  waves  for  ripple-fired  blasts  (e.g.  Gitterman,  1982),  the  main  factors  affecting 
amplitude  of  recorded  seismic  waves  (and  appropriate  magnitude)  are  the  charge  weight 
exploded  simultaneously  (per  delay),  delay  time  and  number  of  delays.  This  analysis  permits 
an  easy  explanation  of  the  observed  differences  for  the  Carlin  and  Kursk  quarry  blasts. 

Fig.  7  plots  magnitude  versus  charge  weight  of  59  ripple-fired  quarry  blasts  in  the 
Galilee,  the  Negev  and  in  Jordan,  showing  considerable  variance  and  no  correlation  for  charges 
of  2-10  tons  and  a  weak  dependence  for  charges  of  10-17  tons.  The  blasting  parameters  used, 
i.e.  delays  of  20-40  msec  between  1 0-40  short,  diagonally  oriented  rows  of  charges  (Gitterman 
and  van  Eck,  1993),  are  close  to  those  at  the  Carlin  mine.  In  this  case  the  blast  size  is  enlarged 
due  to  the  increasing  number  of  delayed  rows  of  hole  charges,  producing  a  prolongation  of 
radiated  energy  and  duration  of  signal  with  no  significant  rise  in  recorded  amplitudes.  It  should 
be  noted  that  local  magnitude  Ml  estimates  based  on  coda  duration  (Shapira,  1988)  are  not 
much  affected  by  this  duration  increase  of  1-1.5  sec.  A  slight  rise  in  magnitude  for  larger  blasts 
is  due  to  some  increase  in  exploded  charge  per  delay. 

Contrary  to  these  observations,  a  clear  correlation  between  charge  weight  and  magnitude 
is  observed  for  10  simultaneous  quarry  shots  (see  Fig.  7)  and  9  Dead  Sea  imderwater 
explosions.  Preliminary  approximate  relationships  are  estimated  (see  Fig.  8): 


for  single  quarry  blasts: 

Ml  =  -2.6+1. 3*log,oW(kg) 

(3.5) 

for  Dead  Sea  UWE: 

Ml  =  -0.01+1  ■2*logjoW(|jg) 

(3.6) 

Data  for  3  Tajikistan  UWE  (kindly  presented  by  Dr.V.Khalturin)  show  similar 
estimates,  as  Dead  Sea  explosions.  A  little  bit  lower  magnitude  values  are  possibly  due  to 
weaker  coupling  (smaller  water  density  and  depth),  and/or  a  different  magnitude  evaluation. 

We  have  also  evaluated  attenuation  of  seismic  waves  as  dependence  of  Peak  Ground 
Velocity  (PGV)  on  the  scaled  distance  R=r(icm/W(Lg)  ,  measuring  maximal  amplitude  on 
seismograms,  for  three  UWE  in  the  Dead  Sea  with  charges  304,  192  and  24  kg,  and  for  10 
single  blasts  (Fig.  9): 


for  single  quarry  blasts: 

PGV,^k„„,  =  4.65.R-'’ 

(3.7) 

for  Dead  Sea  UWE: 

PGV,^t,,e„)  =  67.4.R-'’ 

(3.8) 

The  high  efficiency  of  the  Dead  Sea  UWE  should  be  emphasized  here:  very  small  charges 
of  W=300  kg  yield  high  amplitudes  and  a  local  magnitude  of  Ml=3.1,  which  is  equivalent  to  an 
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Fig.  8.  Seismic  efficiency  of  instantaneous  quarry  blasts  and  underwater  explosions. 
Outstanding  strong  coupling  for  Dead  Sea  UWE  is  observed. 
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Fig.  9.  Empirical  attenuation  of  seismic  waves  for  instantaneous  explosions.  Peak 
Ground  Velocity  versus  scaled  distance. 


instantaneous  quarry  blast  of  more  than  25,000  kg  (Eq.3.5).  According  Eq.  3.4,  a  rough 
estimate  of  a  nuclear  yield  producing  a  similar  magnitude  mb=3.1,  equals  to  16,000  kg.  A 
strong  coupling  (and  minimal  energy  losses),  caused  in  particular  by  the  dense  salt  water  and 
relatively  large  detonation  depth  provides  effective  conditions  for  calibration  explosions 
currently  lacking  in  the  region  and  necessary  for  the  improvement  of  monitoring  the  CTBT. 

The  equations  developed  can  be  used  to  estimate  the  efficiency  of  various  calibration 
explosions  (referring  to  a  signal  amplitude  value  of  0.1  p/sec  as  the  detection  threshold).  For 
example,  for  a  simultaneous  explosion  in  a  quarry,  placed  in  several  dozens  of  spaced  holes,  it 
was  estimated  (Eq.  3.5,3. 7)  magnitude  and  maximal  observational  distance; 

W  =  20,000  kg;  Ml  «  3.0;  *  500  km, 

whereas,  a  smaller  explosion  in  the  Dead  Sea  seems  to  be  much  more  efficient  (Eq.  3. 6,3. 8): 

W  =  5,000 kg:  Ml«4.0-4.2;  «2600km, 

(^max  greatest  distance  at  which  the  event  might  be  observed  on  the  seismograms). 

Similar  UWE  in  the  Dead  Sea  will  probably  be  recorded  by  most  of  the  seismic  stations  and 
local  networks  operating  in  the  Eastern  Mediterranean  and  Middle  East. 

4.  DISCRIMINATION  METHODS  USED  IN  THE  STUDY. 

4.1  Known/Modified  Discriminants. 

4.1.1  Railed  Discriminants:  P/S  ratio.  Coda  wave  statixtics  and  spectral  nvlh 

E/S  ratio  is  the  most  frequently  used  discriminant  in  the  family  of  Pn/Sn,  Pn/Lg  ratios 
which  is  reported  to  be  successful  by  Pomeroy  et  al.  (1982),  Baumgardt  and  Young  (1990)  and 
others.  However,  these  short-period  discriminants  are  heavily  dependent  upon  the  structure  and 
constitution  of  the  uppermost  crust.  Their  performance  varies  firom  region  to  region  and  should 
be  tested  on  a  reliable  database  with  ground  truth  information  in  every  case  study.  Most  of  the 
investigations  mentioned  were  conducted  for  a  single  station  or  beamed  array  data.  Kim  et  al. 
(1994)  showed  that  discrimination  based  on  a  single-station  discriminant  (P/S  spectral  ratio) 
enhances  significantly  when  it  is  averaged  over  the  New  York  State  seismic  network.  We 
obtained  similar  results  in  this  study  by  averaging  amplitude  P/S  and  energy  spectral  ratios 
over  the  ISN. 

It  is  expected  that  the  P/S  ratio  for  explosions  will  be  greater  than  that  for  earthquakes.  In 
practice,  this  appeared  to  be  true  for  only  a  limited  number  of  ISN  stations  at  distances  greater 
than  1 00  km  and  was  realized  only  for  a  small  part  of  the  data  base  with  strong  enough  events. 
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We  tested  this  statistic  in  the  form  of  ratio  of  seismogram  envelope  maximums  Pmax/Smax  in  the 
P  and  S  (Lg,  Rg)  velocity  domains  and  showed  (Gitterman,  et  al.,  1996)  that  only  18  events 
(8  explosions  and  10  earthquakes)  out  of  69  events  in  the  Galilee  with  Ml>2  had  SNR  large 
enough  at  distances  greater  than  100  km.  For  these  data  only  six  remote  sites  out  of  30  stations 
of  the  ISN  had  P/S  values  that  produced  separation  between  earthquakes  and  explosions,  but  at 
relatively  low  efficiency.  Averaging  of  the  ratio  over  these  six  stations  enhances  the  resolving 
power.  No  data  were  available  for  the  remaining  41  events  of  the  data  base,  hence  it  would 
appear  to  be  useless  for  the  most  part  of  our  local  seismicity. 

Coda  wave  mtistics:  Following  the  method  of  Su  etal.(1991),  we  extracted  from  the 
coda  wave  records  (Galilee  dataset)  spectral  amplitude  Ao(f)  and  quality  factor  Q(f) 
characterizing  the  first-order  scattering  coda  wave  model  of  Aki  and  Chouet  (1975).  The 
parameter  Ao(f)  measures  time-independent  spectral  amplitude  of  the  coda,  carrying  spectral 
peculiarities  of  the  source,  van  Eck  (1988)  reported  that  coda  spectra  of  explosions  are  richer  in 
lower  frequency  energy  as  compared  to  earthquakes.  As  shown  in  Gitterman  et  al.  (1996),  the 
ratios  R=Ao(1.5  Hz)/Ao(10  Hz),  averaged  over  ISN  subnetwork  recordings  of  28  earthquakes 
and  30  explosions  in  the  Galilee,  show  much  overlapping  being  only  in  average  higher  for 
explosions  than  for  the  earthquakes.  The  spectral  ratio  of  low-to-high  frequency  energy 
measured  for  the  whole  signal  and  averaged  over  a  subset  of  stations  appeared  to  be  much  more 
efficient. 

Seismograms  of  explosions  usually  have  dominant  surface  waves;  these  waves  propagate 
at  the  upper  crust  layer  with  high  attenuation  owing  to  relatively  low  Q  at  low  frequencies. 
However,  for  the  given  Galilee  dataset  this  method  did  not  provide  separation  for  the  two 
populations,  even  for  the  stronger  events.  Analysis  showed  that  the  method  is  very  sensitive  to 
the  SNR  and,  even  for  regions  where  it  is  applicable,  it  will  be  limited  to  cases  of  SNR  >  10. 

Spectral _ nulls.  Contrary  to  natural  earthquakes,  underground  nuclear  tests  are 

characterized  by  shallow  focal  depths  (several  hundreds  of  meters  and  by  the  spherical 
symmetry  of  the  initial  radiation  from  the  source.  These  features  lead  to  strong  destructive 
interference  between  the  downgoing  P  wave  energy  and  the pP  wave  reflected  from  the  Earth’s 
surface.  This  interference,  however,  is  complicated  by  non-linear  surface  effects  on  the  pP 
reflection  (Lay,  1991). 

The  interference  produces  modulation  minima  in  the  spectra  of  teleseisms  at  frequencies: 

f„,  «  nv/2h  (n  =  1,  2...)  (4.1) 
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where  v  is  the  compressional-wave  velocity  of  the  medium  above  the  source  and  h  is  the  depth. 
For  many  nuclear  tests  conducted  at  depths  h  =  500-700  m  and  for  v  =  1. 1-1.5  km/sec,  Eq.  4.1 
gives  fm  »  1  Hz.  For  deeper  sources  the  velocity  increases  and  the  frequency  value  does  not 
change  a  great  deal.  Spectra  of  the  initial  part  (4-7  sec)  of  the  teleseismic  records  of  nuclear 
explosions  detonated  at  the  Nevada  test  site  demonstrated  pronounced  minima  (nulls)  near 
1  Hz,  which  could  be  interpreted  in  terms  of  interference  of  P  and  />P-wave  (Kulhanek,  1971). 

Being  a  source-like  effect,  the  spectral  null  is  apparent  at  about  the  same  frequency  in 
spectra  of  different  stations  of  a  network  and  can  be  used  to  discriminate  nuclear  tests  from 
deep  earthquakes.  This  phenomenon  was  also  observed  for  ISN  short-period  (Fig.  10)  and 
broad-band  (Fig.  11)  records  of  several  Chinese  tests,  with  almost  identical  deep  spectral 
minima,  which  suggests  that  the  locations  and  depths  of  these  explosions  are  very  close.  (An 
identity  of  SP  and  BB  seismograms  and  spectra  should  be  noted,  thus  confirming  absence  of 
significant  distortions  in  low-frequencies  of  teleseismic  P-waves  recorded  by  short-period 
instruments). 

We  tried  to  check  the  azimuthal  independence  inherent  to  this  source-like  effect  and 
compare  spectra  of  the  P-waves  group  at  different  ISN  stations.  Some  explosions  (e.g.  NT  15 
on  Fig.  12)  demonstrate  a  distinct  spectral  modulation  and  high  coherency  of  spectral  shapes  at 
many  stations,  with  minima  at  about  1  Hz  (several  multiple  minima  can  be  observed  at  2,  3,  4 
and  5  Hz).  Nevertheless  many  explosions  revealed  strong  variability  of  the  null  frequency 
(Fig.  13),  which  can  hardly  be  identified  for  some  stations.  On  the  other  hand,  several 
earthquakes  (e.g.  QT26)  also  show  the  same  spectral  null  phenomenon,  possibly  caused  by 
interference  of  rupture  sub-events  (see  Fig.  14).  We  conclude  that  this  spectral  feature  cannot 
be  utilized  here  as  a  reliable  discriminant. 


Successful  Discriminants:  ener^  spectral  ratio,  spectmeram  and  P/P-cnda  rntin 
Energy  spectral  ratio  (SR).  ISN  seismograms  of  quarry  blasts  demonstrate  more  low- 
frequency  content  of  seismic  waves  than  local  earthquakes  (Fig.  15).  It  has  been  suggested  that 
spall  mechanisms  contribute  to  low  frequency  seismic  energy  and  the  general  tendency  for 
quarry  blasts  to  appear  deficient  in  high  frequency  energy  (Barker  et  al.,  1993).  We  suggest  that 
the  phenomenon  is  mainly  due  to  dominant  surface  waves  (including  Rg  on  the  vertical 
component)  generated  by  quarry  blasts  and  associated  with  the  regional  crustal  structure  of 
widespread  unconsolidated,  low  velocity  and  low  Q  subsurface  sediments  (e.g.  Ginzburg  and 
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Fig.  10.  Several  Chinese  tests  recorded  at  the  ISN  station  PRNI  (a).  Instrument-corrected 
spectra  of  time  windows  about  7  sec,  marked  by  arrows,  show  distinct  spectral  nulls 
(b),  almost  identical  for  all  events.  After  enhancement  of  spectra  resolution  by 
expanding  the  analyzed  signal  by  zeros  to  28  sec,  and  reducing  Af  from  0.1  Hz  to  0.02 
Hz,  the  spectral  nulls  became  deeper  and  better  expressed  (c). 
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FFT  from  datasets  :  1 ,  2,  3 ,  4 •  Sm  Wd  :  I 


Fig.  11.  Records  and  spectra  of  the  Chinese  tests  at  the  BB  station  BGIO  (vertical 
component)  demonstrating  similarity  to  ISN  short-period  data  The  analyzed 
signal  m  time  windows  about  7  sec,  marked  by  arrows,  was  expanded  by  zeros  to 
28  sec,  for  better  spectra  resolution. 
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Nuclear  test  NT15  88/08/22,  W.Siberia,  66“N  78.5°E  ,  mb=5.3 


FFT  spectra  at  ISN  stations,  time  window  [0-15  sec] 


Fig.  12.  The  nuclear  test  NT15,  demonstrating  distinct  spectral  nulls  and  modulation  and 
high  coherency  of  spectral  shapes  at  different  stations. 
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Nuclear  test  NT27  95/06/15,  China  (Lop  Nor),  42®N  89®E,  mb=6-1 
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Fig.  13.  The  nuclear  explosion  NT27  showing  strong  variability  of  the  null  frequency  and 
spectral  shapes.  ^ 


Earthquake  QT26  91/08/19,  China,  47®N  85°E,  H=44  km,  mb=5.7 


SRGI  Rug  19.91  (06:12)  - - - - - — - - - — _ _ 

(1.155)  (1) 

(1.1S4)  Cl) 

lyM/A 

(1.000)  (I) 

va|Av/ 

•  . 

j\  (1.0Sl>  (1) 

.11  ^  3.2  34  36  38  40  43  44  46  48  S0  53  ^  ^=) 

FFT  spectra  at  ISN  stations,  time  window  [0-15  sec] 


Fig. 


14.  The  earthquake  QT26,  revealing  explosion-like  nulls  about  1  Hz,  modulation  and 
hgh  coherency  m  spectra  for  19  stations.  Seismograms  of  four  stations  only  are 
presented  for  better  observation  of  characteristic  interference  features  in  the  P-wave. 


Fig.  15.  Examples  of  recordings  at  different  ISN  stations  and  appropriate  spectra:  (a)  a 
ripple-fired  blast  (ES6)  fi'om  a  Negev  quarry,  showing  distinct  azimuth  and 
distance-independent  spectral  modulation,  and  low-frequency  spectral  content 
due  to  strong  surface  (Rg)  waves. 
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Fig  15b.  An  earthquake  from  the  Dead  Sea  (QS3),  demonstrating  non-coherency  of 
spectral  shapes  and  high-frequency  energy.  Arrows  show  analyzed  time 
windows. 
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Folkman,  1980),  trapping  low  frequency  energy  and  absorbing  high  frequencies  along 
propagation  paths.  This  explanation  is  supported  by  similar  observations  from  New  England 
(Li  etal.,  1996)  and  by  the  reverse  situation  for  quarry  blasts  in  the  Scandinavian  shield 
(Baumgardt  and  Young,  1990). 

We  utilized  this  effect  in  the  energy  spectral  ratio  of  seismic  energy  in  the  low- 
frequency  range  [f,,f2]  and  the  high-frequency  range  [f3,f4]  (Gitterman  and  Shapira,  1993): 

=  I  P(/)prf//j  |5(/)|V/  (4  2) 

/.  h 

where  |S(f)|  is  the  smoothed  spectrum  of  ground  velocity  for  the  whole  signal.  This  ratio, 
implying  evident  physical  meaning,  is  similar  to  the  ratios  of  average  spectral  amplitude 
(Beimett  and  Murphy,  1986,  Walter  etal.,  1995),  or  power  spectrum  (Pulli,  1995),  calculated 
for  specific  regional  phases.  For  an  event,  values  at  different  stations  show  strong 
scattering;  averaging  over  a  subnet  is  crucial  for  separation  between  earthquakes  and 
explosions  (Shapira  et  al.,  1996). 

Spectrogram  is  a  robust  characterization  of  the  seismic  energy  distribution  across  the 
specified  frequency  interval,  often  useful  in  the  multivariate  discrimination  analysis  (Dowla, 
1995,  Pulli,  1995).  It  is  a  vector  statistic  resulting  from  calculation  of  RMS  of  FFT  amplitudes 
in  the  series  of  spectral  windows.  For  the  local  data  analysis  spectrograms  were  calculated  at 
the  large  enough  time  intervals  comprising  together  all  seismic  phases:  P,  S,  Lg,  Rg  and  a  part 
of  coda.  For  teleseismic  discrimination  only,  P  and  P-coda  time  intervals  were  considered. 

P/P-coda  ratio  teleseismic  discriminant  resembles  classic  complexity  statistics  (Dahlman 
and  Israelson,  1977)  and  was  shown  to  be  efficient  for  the  NORESS  teleseismic  recordings 
(Tsvang  et  al.,  1993). 

P/P-coda  ratio  results  from  the  following  computations: 

•  Spectra  are  computed  for  the  forward  short  time  interval  (2  seconds)  and  the  following 
long  time  interval  (30  seconds). 

•  A  spectrogram  is  computed  for  each  of  the  time  intervals  in  several  (six)  equal  frequency 
bands  in  specified  spectral  interval  (from  0  to  5  Hz). 

•  The  maximum  coefficient  is  taken  in  each  of  the  two  spectrograms  and  averaged  over  a 
set  of  the  network  stations.  The  resulting  statistics  are  used  to  calculate  the  P/P-coda  ratio 
discriminant. 
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4.2  Newly  Developed  Procedures. 

4-2.1  Snectral  Semblance  (SS).  Distinctive  features  of  the  quarry  blasts  and  underwater 
explosions  records  are  clear  spectral  modulation  patterns  presented  by  coherent  minimfl  (or 
nulls)  and  maxima  observed  at  low  frequencies  1-10  Hz  (Fig.  15a).  This  azimuth-independent 
modulation  is  caused  by  constructive-destructive  interference  of  linearly  superimposed 
waveforms  from  spaced  charges  exploded  at  millisecond  delays  (“ripple-firing”)  (e.g. 
Gitterman,  1982;  Stump  and  Reinke,  1988),  or  from  oscillations  of  gaseous  bubbles  in  water 
and  indicative  of  the  source-effect.  Another  approach  to  reveal  spectral  modulation  involves  the 
time-independent  pattern  in  a  sonogram  at  a  single  station  (e.g.  Hedlin  et  al.,  1989),  but  this 
kind  of  modulation  can  appear  also  due  to  a  propagation  path  effect. 

In  contrast  to  a  ripple-fired  blast,  the  spectral  pattern  for  an  earthquake  is  strongly 
azimuth-dependent,  owing  partly  to  the  directivity  effect  (e.g.  Bakun  et  al,  1978),  therefore,  for 
stations  with  broad  azimuthal  coverage,  we  observe  an  irregular  character  of  spectral  shapes 
and  minima  for  different  azimuths  (Fig.  15b). 

The  semblance  and  cross-correlation  statistics  are  commonly  used  in  seismic  prospecting 
for  phase  correlation  of  seismic  traces  in  the  time  domain  (Neidell  and  Taner,  1971).  We 
modified  these  statistics  to  assess  quantitatively  the  coherency  of  spectral  shapes  for  different 
stations.  As  expected,  a  clear  equivalence  between  the  two  characteristics  was  observed 
(Gitterman  etal.,  1996),  therefore,  in  the  following  we  considered  only  spectral  semblance 
(SS),  i.e.: 
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where  S^;  =  logioSk(fj)  -  log  spectral  amplitude  at  k  station;  S,j  is  the  average  spectral  level;  N  is 
the  number  of  used  stations  and  [Fi,F2]  is  the  frequency  band  for  the  analysis. 

4.2.2  Velogram  Analysis  Statistics.  Another  new  RDSN  oriented  discriminant  is  based  on 
velograms  -  signal  envelopes  versus  group  velocity.  The  signal  envelopes  are  calculated  using 
RMS  in  the  moving  window,  and  group  velocities  are  the  ratios  of  R/T,  where  R  is  the  distance 
and  T  is  the  lapse  time.  The  analysis  is  based  on  computation  of  the  group  velocity  value 
corresponding  to  the  local  maximum  of  the  velogram  in  the  range  of  1-4  km/sec  (see  Fig.  16). 
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The  relationship  between  the  S-waves  group  velocity  and  the  distance  R  (in  the  range 
10-150  km)  differs  for  blasts  and  earthquakes  (Fig.  17).  The  effect  should  probably  be 
attributed  to  the  different  excitation  of  surface  waves  from  these  two  types  of  seismic  events 
or/and  different  S-waves  group  velocities  for  shallow  and  deep  sources  emphasized  in  many 
studies  (Kafka,  1990;  Walter  etal.,  1996;  Alexander  etal.,  1995).  It  was  shown  (Pinsky  and 
Shapira,  1997),  that  a  simple  “C’  statistic,  derived  from  the  parametric  approximation 

Vms(R)  =  a  +  b*ln(R)  (4.4) 

C=b+0.33*a,  (4.5) 

provides  the  true  identification  of  the  events  from  the  Galilee  data  base  recorded  by  the  ISN 
stations. 

4.3  Multivariate  Procedures  of  Pattern  Recognition 

An  additional  effect  is  achieved  using  an  Integrative  Approach  (lA)  using  a  number  of 
different  physically  signal  features.  Wuster  (1993)  showed  that  even  the  simple  majority  voting 
of  several  different  physical  discriminants  with  equal  weights  provides  improvement  of 
classification  result. 

4.3.1  Linear  Discrimination  Function  (LDF)  -  commonly  used  in  discrimination  seismology 
statistical  procedure  which  accounts  for  different  weights  of  discriminants  and  their  mutual 
correlation.  It  is  based  on  the  assumption  that  the  vectors  X  of  observed  parameters  generated 
by  two  stochastic  (physical)  mechanisms  H,,  H2  (earthquakes  and  explosions),  are  described  by 
Gaussian  distribution  with  equal  covariance  matrices  S=Si  =82  and  different  mathematical 
expectations  Mj.  The  optimal  decision  rule,  based  on  log-likelihood  ratio  function  V(X),  is 
(Tsvang  et  al.,  1993): 

(assign  X  to  H,,  if  V(X)  >  0  and  assign  X  to  H2,  if  V(X)  <  0}  (4.6) 

where  V(X)  =  F,(X)  -  F2(X)  =  2(M2  -  M^^S'^X  -  M2'^S''M2  is  a  linear  function  of 

vector  X. 

Practically,  instead  of  covariance  matrix  Sj  and  expectation  M;  corresponding  statistics 
are  substituted  based  on  appropriate  data  samples.  The  assumption  that  Sj  =  82  ,  yielding 
separation  of  the  two  data  sets  by  the  hyperplane  V(x)=0,  is  rarely  realized  and  the  general 
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Fig.  17.  Velogram  analysis:  Vms-distance  ciuves  for  an  earthquake  and  an  explosion. 
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expression  for  the  likelihood  function  should  be  used  for  optimal  decisions.  However,  in 
practice  even  in  this  case,  the  algorithm  usually  works  fine  if  sample  covariances  Sj  and  S2  are 
substituted  by  their  average  value. 

4.3.2  Artificial  Neural  Network  fANN^  should  have  an  advantage  over  the  LDF  in  the  case 
where  the  separating  surface  between  the  data  sets  is  essentially  nonlinear. 

We  have  used  the  ANN  program,  which  was  kindly  provided  by  Dr.  F.  Dowla  (LLNL, 
USA).  This  network  is  a  classical  example  (Dowla,  1995)  of  a  supervised  learning  network 
applying  backpropagation  learning  algorithm  to  associate  inputs  with  corresponding  outputs  for 
all  or  most  of  the  events  in  a  training  set.  This  algorithm  uses  a  gradient  descent  method  to 
systematically  modify  the  weights  in  the  network  so  as  to  minimize  the  network  output  error. 
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5.  ANALYSIS  OF  DISCRIMINANT  PERFORMANCE 


5.1  Application  of  Energy  Spectral  Ratio  and  Spectral  Semblance  to  Regional  and 

Teleseismic  ISN  and  NORESS  Recordings. 

Regional  data.  Initially,  using  data  from  the  Galilee  region,  we  tried  for  spectra  radio 
fi=l  Hz,  f2=f3=6  Hz,  f4=ll  Hz  and  Fj  =1  Hz,  F2=12Hz  for  spectral  semblance  (Gitterman 
etal.,  1996),  thus  facilitating  almost  the  whole  recording  frequency  band.  Application  of 
another  ratio  (2-4 Hz)/(6-8  Hz),  used  by  Li  etal.  (1996)  for  the  Lg  spectra  to  separate  small 
events  in  New  England,  reduced  the  resolving  power  (Fig.  18a).  However,  after  a  low 
frequency  shift,  the  (1-3  Hz)/(6-8  Hz)  ratio  enhanced  separation  between  quarry  blasts  and 
earthquakes  and  provided  the  best  discriminant  performance  (Fig.  18b).  It  is  interesting  to  note 
that  almost  the  same  frequency  bands  (1-2  Hz)/(6-8  Hz),  but  calculated  for  separate  seismic 
phases,  were  used  for  discrimination  of  NTS  explosions  and  earthquakes  at  short  distances,  and 
performed  well  over  magnitude  3.5  (Walter  et  al.,  1995). 

As  expected,  the  semblance  values  are  much  higher  for  quarry  blasts  than  for  earthquakes 
and  remain  practically  unchanged  when  using  other  (different  from  usual  0.5  Hz)  windows  of 
spectra  smoothing  (0.25  and  1  Hz). 

For  the  few  earthquakes  originating  on  the  main  fault  in  the  Sea  of  Galilee,  we  obtained 
anomalous  high  SR  and  SS  values.  A  different  range  (1-7  Hz)  provided  “normal”  semblance 
values  for  these  earthquakes,  but  for  the  rest  of  events  we  obtained  a  worse  separation.  The  best 
semblance  performance  for  the  Galilee  and  Gilad  regions  was  obtained  for  the  (1-8  Hz)  band 
(Fig.  18c),  nevertheless,  due  to  worse  results  for  quarry  blasts  and  earthquakes  from  the 
southern  Dead  Sea/Negev  area  and  single  blasts  (Fig.  1 8d),  the  optimal  range  remained  valid  as 
1-12  Hz.  In  general,  all  earthquakes  show  consistent  ratio  and  semblance  values,  lower  than 
those  for  ripple-fired  quarry  blasts  with  the  thresholds:  Re  =6.1,  Sf  =  0.69. 

The  discrimination  results  for  all  local  events  are  presented  in  Fig.  19.  Anomalous  low 
ratio  values  for  Dead  Sea  UWE  are  possibly  associated  with  the  higher  frequency  content  of 
radiated  seismic  waves  due  to  deeper  source  (70m)  than  that  of  the  off-shore  UWE  near  Tyre 
(assumed  10-30m),  higher  density  of  water  and  small  charges  (mostly  16  and  24  kg). 
According  to  Eq.  3.1,  for  the  charge  weight  W=24  kg,  the  predicted  dominant  frequency  due  to 
the  bubble-effect  is  fib«6.7  Hz  and  is  close  to  the  observed  value  of  6.0-6. 5  Hz  (see  Fig.  20). 
This  observation  could  also  be  influenced  by  the  fj  harmonic  (6.3  Hz  from  Eq.  3)  related  to  the 
detonation  depth.  It  turned  out  that  for  small  explosions  the  reverberation  effect  caused  by  the 


Spectral  semblance  (1-12  Hz) 


Fig.  18.  Discrimination  results  (spectral  semblance  vs  spectral  ratio)  for  local  events  at 
different  frequency  ranges.  Separation  of  earthquakes  and  quarry  blasts  in 
Galilee  was  improved  after  a  shift  of  low-frequency  band  in  the  energy  ratio 
from  (2-4Hz)  (a)  to  (l-3Hz)  (b);  the  semblance  range  (l-8Hz)  provided  the 
best  separation  of  the  events  from  Galilee  and  Gilad  regions  (c),  but  showed  a 
reduced  resolving  power  for  the  Southern  dataset  and  single  quarry  blasts  (d). 
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Spectral  semblance  (1-8  Hz) 


Fig.  19.  Discrimination  results  (spectral  semblance  vs  spectral  ratio)  for  all  local  events 
at  the  optimal  frequency  ranges  (a).  The  semblance  frequency  range  (1-8  Hz) 
produces  worser  separation  and  stronger  overlapping  of  earthquaJkes  and 
explosions  (b). 
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water  depth  H  (and  the  appropriate  frequency  y  was  not  revealed  and,  therefore,  did  not 
contribute  to  the  low-frequency  energy  in  the  band  1-3  Hz  (unlike  the  largest  EXl,  see  Fig.  5). 
We  also  suggest  that  compressional  waves  from  small  explosions  are  scattered  stronger  (and 
lose  energy)  during  reverberational  multiple  reflections  from  the  uneven  bottom,  due  to  smaller 
wave  lengths  (A,=Vo/fb=l 770/6.7=260  m),  as  compared  to  the  explosion  EXl  (;\,=610  m). 

The  spectral  ratios  for  seven  instantaneous  (single)  blasts  in  the  Mt.  Nitzim  quarry  are 
low  (R£=2-5),  unlike  those  for  three  single  shot  blasts  at  the  Revaya  quarry  (Re>100), 
demonstrating  strong  low-frequency  seismic  radiation  (Fig.  21).  This  is  possibly  associated 
with  the  uppermost  geology  near  the  two  quarries  located  in  different  regions  (see  map  in 
Fig.  1).  The  Revaya  quarry  is  located  in  the  Beit-Shean  sedimentary  basin  formed  by  soft  river 
and  lake  deposits,  and  the  propagation  paths  to  the  stations  cross  the  basin,  producing  the  low 
frequency  surface  waves  due  to  channeling  and  the  intrinsic  absorption  of  high  frequency 
energy.  In  contrast,  the  Mt.  Nitzim  site,  containing  consolidated  massive  carbonates,  is  more 
correspondent  to  the  hardrock  conditions  producing  high  frequency  seismic  radiation  (Fig.22). 

The  semblance  values  for  single  explosions  are  lower  than  those  for  ripple-fired  blasts, 
but  higher  than  for  earthquakes  (see  Fig.  19).  A  possible  explanation  for  this  relates  to  the  fact 
that  single  shots  do  not  possess  the  azimuth-independent  modulation  in  spectra  caused  by  time 
delays  and  producing  high  semblance  values.  Nevertheless,  they  are,  actually,  point  sources, 
and  dependence  of  radiated  waves  on  azimuth  is  less  significant,  than  for  earthquake  line 
sources  subjected  to  the  directivity  effect.  Spectra  at  ISN  stations  from  single  explosions  show 
some  coherent  gross-structure  (envelope  of  random  small  scallops):  a  broad  maximum  at  low 
frequencies  (1-5  Hz),  mainly  due  to  surface  waves,  and  steep  descend  at  high  frequencies 
(5-12  Hz),  caused  by  intrinsic  attenuation  (see  Fig.  21,22).  This  gross-structure  provides  a 
significant  contribution  to  higher  (than  for  earthquakes)  semblance  values.  This  effect  is 
inherent  also  in  spectra  of  ripple-fired  blasts. 

Tclcssisniic — dp. to..  The  same  SR  and  SS  procedures  were  applied  to  68  Eurasian 
teleseismic  events,  with  the  only  evident  correction  (compared  with  local  events)  being  the 
shift  of  spectral  bands  to  lower  frequencies. 

The  analyzed  waves  {P  and  P-coda)  from  teleseismic  earthquakes  have  more  low- 
frequency  energy  and  more  coherent  spectral  shapes  than  nuclear  explosions  (see  Figs.  12-14). 
Therefore,  the  semblance  and  spectral  ratio  values  of  ISN  records  of  the  earthquakes  are  higher 
than  those  for  the  explosions  (see  Fig.  23),  which  is  opposite  to  observations  for  local  events 
(see  Figs.  18  and  19).  There  is  a  number  of  short-period  observations  of  Soviet  and 
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Fig.  21.  Three  single  blasts  at  Revaya  quarry  (SB1,SB2,SB3,  W=1946-2840kg,  r=3- 
29  km),  with  prominent  surface  waves,  showing  a  high  identity  of  wave  forms 
(a)  and  spectral  shapes  (b)  at  station  GVMR  (analyzed  time  window  t «  20  sec). 
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Fig.  22.  Two  single  blasts  at  Mt.Nitzim  quarry  (SB4,  W=2400kg,  SB5,  2200kg,  r=22- 
35kra)  demonstrating  non-regular  wave  forms  (a)  and  hi^-frequency  spectra 
(b)  at  different  stations  (analyzed  time  window  t «  20  sec). 

40 


Eurasian  teleseismic  events  (ISN  records,  normins) 


29  nuclear  tests  fc),  39  earthquakes 


Spectral  semblance  (0.6-2Hz) 
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Fig.  23.  Discrimination  results  (spectral  semblance  vs  spectral  ratio,  at  optimal 
frequency  ranges)  for  teleseismic  Eurasian  earthquakes  and  nuclear 
explosions,  based  on  ISN  recordings,  without  (a)  and  with  (b)  remove 
instrument  response. 
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Chinese  nuclear  tests  at  different  sites,  confirming  more  high-frequency  energy  for  the 
explosions  than  for  the  shallow  crustal  earthquakes  (e.g.  Tailor  and  Marshall,  1991,  for  UK 
teleseismic  arrays).  This  can  be  explained  by  a  smaller  spatial  area  of  energy  release  for  nuclear 
tests  (the  finite  elastic  radius)  and  shorter  duration  of  the  source  function  (the  time  rise),  as 
compared  to  earthquakes  (Savage,  1972),  which  produces  a  shorter  wave  length  of  P-waves. 
Azimuths  of  a  teleseismic  event  to  the  stations  are  very  close  (in  contrast  to  local  earthquakes), 
thus  significantly  increasing  spectra  coherency.  Nevertheless,  short-period  explosion  signals 
(including  scattered  waves  in  P-coda),  are  more  influenced  by  inhomogeneities  in  the  upper 
crust  beneath  the  ISN  stations  spaced  several  dozen  kilometers  apart,  producing  a  decreasing  of 
spectra  coherency. 

The  best  discrimination  performance  was  provided  in  the  frequency  range 
(0.6-1  Hz)/(l-3  Hz)  for  averaged  spectral  ratios  and  in  (0.6-2  Hz)  for  semblance,  calculated 
over  10-15  ISN  stations.  It  turned  out  that  removing  instrument  response  from  the  spectra 
significantly  improves  the  spectral  semblance  performance  (Fig.  23  a,b).  Joint  application  of 
the  two  discriminants  provided  good  separation  between  the  earthquake  and  explosion 
populations  with  a  few  events  overlapping. 

We  tried  to  track  a  possible  dependence  of  SS  and  SR  statistics  on  magnitude  and  depth 
(for  earthquakes  only).  The  results  are  presented  on  Figs.  24  and  25.  The  only  relationship 
observed  is  between  the  spectral  ratio  and  the  magnitude  (Fig.  24a),  which  can  be  explained  by 
increasing  the  wavelength  and  low-ffequency  component  of  radiated  energy  for  larger  events. 
A  visible  saturation  of  this  dependence  for  mb>6  is  possibly  due  to  a  shift  of  a  fraction  of  low- 
fi-equency  energy  beyond  the  high  pass  cut  off  (0.2Hz),  or  the  analysis  band  (0.6-3Hz). 

In  the  NORESS  data  we  foimd  a  similar  high-frequency  character  of  the  P-arrivals  for 
explosions  relative  to  earthquakes  (see  Figures  26,27),  consistent  with  results  of  Tsvang  et  al. 
(1993).  Therefore,  spectral  ratios,  which  were  calculated  (without  removing  instrument 
response)  for  8-10  stations  placed  over  the  entire  array  (Fig.  26c),  showed  the  same  (as  for 
ISN)  interrelation  between  the  two  populations,  whereas  the  character  of  the  semblance 
changed  to  the  opposite,  i.e.  very  high  (Sf  >0.94)  close  values  for  explosions  and  lower,  in¬ 
average,  dispersed  values  for  earthquakes  (Fig.  28).  In  this  case,  owing  to  the  small  dimensions 
of  the  array  (diameter  about  3  km),  local  inhomogeneities  in  the  upper  crust  and  azimuth 
differences  are  negligible.  We  therefore  suggest  that  the  main  factor  affecting  station  spectra 
coherency  is  a  waveform  character  in  the  0-15  sec  window:  a  simple  predominant  P-signal 
pulse  for  explosions,  and  several  equal-amplitude  signals  of  different  shape  for  earthquakes. 


42 


EURASIAN  TELESEISMIC  EVENTS  (ISN  records) 


10  :t 


100  150 

Depth  (km) 


EURASIAN  TELESEISMIC  EVENTS  (ISN  records) 


100  150 

Depth  (km) 


Fig.  25.  Spectral  ratio  (a)  and  spectral  semblance  (b)  vs  depth  for  ISN  records  of 
teleseismic  earthquakes. 
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Fig.  26.  Typical  seismograms  (a)  and  spectra  (b)  of  a  nuclear  test  (NTH)  recorded  at 
NORESS.  The  stations  used  for  the  analysis  are  placed  on  all  array  rings  (c). 
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28.  Discrimination  results  (spectral  semblance  vs  spectral  ratio)  for  teleseismic 
earthquakes  and  nuclear  tests,  recorded  by  NORESS.  Different  frequency  bands  for 
spectral  ratio  are  used:  (0.6-0.8Hz)/(l-3Hz)  (a),  and  (0.6-0.8Hz)/(0.6-10Hz)  (b). 
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Joint  application  of  the  two  statistics  with  the  same  as  for  ISN  frequency  bands 
(0.6-1  Hz)/(l-3  Hz)  provided  reliable  discrimination  between  the  two  types  of  seismic  events, 
with  the  only  outlier  being  the  deep  Afghanistan  earthquake  QUIO  (Fig.  28a).  By  applying 
different  frequency  bands  (0.6-0.8  Hz)/(0.6-10  Hz),  which  is  actually  the  normalized  low- 
frequency  energy  in  the  very  narrow  band,  we  obtained  full  separation  of  the  earthquake  and 
explosion  populations  (Fig.  28b). 

In  general,  the  spectral  ratio  performance  appears  to  be  better  for  NORESS  recordings, 
whereas  the  semblance  is  better  for  ISN  data. 

5.2  Automated  classifiers  for  local  and  teleseismic  ISN  records. 

Automated  Signal  Classification  (ASC)  is  a  significant  challenge  in  the  monitoring  of  a 
CTBT.  We  constructed  a  reliable  automated  classifier  for  regional  and  teleseismic  events 
identification  as  a  combination  of  the  most  promising  discriminants  and  using  the  advantage 
provided  by  the  ISN  as  a  regional  dense  seismic  network  (RDSN). 

Method.  We  select  several  useful  parameters,  choose  empirical  thresholds  and  evaluate 
rate  of  errors,  integrated  by  the  majority  voting  scheme.  Simple  majority  voting  appears  to  be 
useful  for  multidimensional  classification  analysis  of  local  events,  but  it  does  not,  however, 
account  for  the  mutual  dependency  (cross-correlation)  and  different  weights  of  scalar  statistics 
and  may  be  regarded  only  as  the  first  order  approximation  of  error  probability.  We  therefore 
address  more  sophisticated  classifiers  such  as  the  LDF  and  ANN,  applying  them  to  some  of 
the  most  promising  combinations  of  parameters  and  evaluate  the  error  probability  through  the 
leave-one-out  cross-validation  test  (Fukunaga  and  Hummels,  1989)  to  create  an  independent 
test  data  base.  The  flow  chart  of  the  multivariate  classification  used  in  the  study  is  presented  on 
Figure  29. 

Local  training  data  set  for  construction  of  ASC  comprises  80  quarry  blasts  and 
65  earthquakes  with  a  magnitude  range  Ml=1  .5-2.5  and  a  depth  range  H=0-28  km,  from  three 
geologically  different  non-overlapped  regions  of  Israel  and  Jordan:  the  Galilee,  southern  Dead 
Sea  and  Gilad  (see  map  in  Fig.  1  and  Table  1).  These  data  are  part  of  the  general  data  base 
described  in  Chapter  2.1.  The  events  are  characterized  by  enhanced  excitation  of  surface  waves 
for  explosions  as  compared  to  earthquakes.  This  is  possibly  the  main  reason  for  the 
predominant  low  frequency  and  coherency  of  spectral  shapes  at  different  stations,  as  well  as  the 
lagging  behind  of  the  S-wave  trains  for  explosions  as  compared  to  earthquakes. 

These  three  phenomena  have  been  quantified  using  the  following  discriminants: 
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Fig.  29  Flow  chart  of  the  multivariate  classification 
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•  spectrogram  for  different  partition  of  the  (1-11  Hz)  frequency  band.  Two  types  of 
spectrogram  were  considered,  depending  on  the  normalization  coefficient: 
“spectrogram  SI ’’-normalization  by  the  maximum  of  the  spectrogram  coefficients; 

spectrogram  S2  -  by  the  squared  vector  norm  (signal  energy).  Separate  coefficients 
of  the  spectrogram  vectors  were  also  considered. 

•  energy  spectral  ratio  -  low-/high-frequency  ratio,  averaged  over  a  set  of  stations; 

•  spectral  semblance  -  coherency  measure  of  spectral  shapes  for  a  set  of  stations; 

•  velogram  parameter  C  -  kinematic  scalar  discriminant. 

All  the  discriminants  mentioned  integratively  utilize  the  multistation  character  of  the 
recordings.  Combinations  of  the  discriminants  are  presented  in  Table  3.  Thresholds  for 
individual  statistics,  separating  blasts  and  earthquakes,  were  chosen  empirically,  to  minimize 
the  number  of  errors.  The  majority  voting,  LDF  and  ANN  procedures  were  applied  to  the 
different  combinations  of  statistics.  The  number  of  errors  for  LDF  and  ANN  was  estimated  by 
the  leave-one-out  test. 


Table  3.  Discrimination  results  for  the  local  training  dataset  by  Multivariate  Procedures 


Parameter 

^Statistic/  spectral 
Band 

Threshold 

Level 

Number  of  Errors 

Voting 

leaveHMie^ouMest 

LDF 

ANN  ■ 

N1 

log,  spectral  ratio, 

(1-3  Hz)/6-8  Hz) 

1.8 

2 

2 

3 

N2 

spectral  semblance, 
(1-12  Hz) 

0.7 

4 

8 

6 

N3 

velogram  parameter, 

C  (1-10  Hz) 

0.68 

4 

9 

8 

N4 

spectrogram  SI  coef 
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Results  of  the  local  events  classification  using  the  four  scalar  statistics  specified  in 
Table  3  are  presented  in  Fig.  30  together  with  depth  specification  at  the  top  and  the  results  of 
majority  voting  at  the  bottom.  The  data  in  Fig.  3  are  split  according  to  regions  of  the  event 
locations.  The  Voting  is  defined  as  the  sum  of  imits  for  the  “correct  vote”  and  zeroes  for  the 
“wrong  vote”;  hence  a  “correct”  majority  decision  corresponds  to  the  Voting>2  for  the  four 
voters  case.  The  results  of  the  analysis  of  Fig.  30  are  as  follows: 

•  The  statistics  selected  exhibit  from  good  (four  errors  for  the  semblance  and  the  velogram 
parameter  Q  to  excellent  (two  errors  for  the  spectrogram  7-9  Hz  and  the  spectral  ratio)-, 
144  events  out  of  145  are  correctly  classified  using  the  majority  Voting  >2  criterion. 

•  The  worst  Voting=\  is  obtained  for  event  N1 19  (QS9)  -  an  earthquake  from  the  Dead  Sea 
region  with  zero  depth  specification.  The  only  positive  voice  for  it  belongs  to  the  velogram 
parameter  C. 

•  1 1  events  have  Voting=3  value  and  are  correctly  classified.  Only  one  of  these  earthquakes, 
event  N120  (QSIO),  may  be  specified  as  shallow. 

Fig.  31  A,B  presents  the  results  of  the  Error  probability  estimation  by  the  leave-one-out 
test  for  the  LDF  and  ANN,  applied  to  the  statistics  described  in  Table  3.  Fig.  3 1  shows  that  the 
error  probability  for  the  classifiers  varies  between  0.6%  to  6.2%,  where,  on  average,  ANN 
performs  better  than  LDF.  The  best  computational  scheme  is  ANN  +  spectrogram  N1  (5  band 
case),  which  provided  only  one  mistake  for  the  Dead  Sea  earthquake  (QS9).  The  second  place 
with  an  error  probability  of  1 .4  %  shares  the  ANN  +  spectrogram  N1  (7-9  Hz)  (relative  energy 
in  the  7-9  Hz  band)  and  LDF  +  Energy  Spectral  ratio  -  both  scalar  statistics. 

Teleseismic  training  data  set  is  described  in  Chapter  2.2.  The  seismograms  of  the 
nuclear  test  are  characterized  by  relatively  higher  frequencies  of  particle  motion  and  impulse 
shape  of  the  P-wave  onset.  Most  of  the  explosions  (18)  are  from  the  Shagan  river 
(Semipalatinsk)  test  site,  eight  from  China  (Lop  Nor),  one  from  Western  Siberia  and  two  from 
Novaya  Zemlya  test  sites.  The  nuclear  detonations  took  place  in  various  geological  conditions 
and  show  significantly  different  waveforms,  sometimes  similar  to  those  of  earthquakes.  The 
deep  events,  mostly  from  northern  Afghanistan,  exhibit  usually  impulsive  P  wave  seismograms 
similar  to  most  of  nuclear  tests  mentioned  above  (see  Fig.  32). 

On  this  data  set  we  investigated  the  combinations  of  P/P-coda  ratio,  spectral  semblance 
(0.6-2  Hz)  and  energy  spectral  ratio  (0.6-1  Hz)/(  1-3  Hz)  statistics.  The  last  two  spectral 
statistics  were  investigated  in  two  modes:  with  and  without  prefiltering  of  the  seismograms  by 
the  (0.5-5  Hz)  Butterworth  filter. 
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Fig.  30.  Majority  voting  for  the  four  main  scalar  statistics  (parameter  numbers  from 

^ble  3),  applied  to  the  ISN  records  of  80  quarry  blasts  (0)  and  65  earthquakes 
(+)  from  Gahlee,  Dead  Sea  and  Gilad  region  (local  training  data  set). 
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Parameter  number 


ANN  Performance 


Parameter  number 


Parameter  Parameter 


Fig.  31.  Error  probability  estimation  by  the  leave-one-out  test  of  the  LDF  and  ANN 
procedures,  applied  to:  (A,B)-local  training  data  set  (parameter  numbers  from  Table 
3);  (C,D)-teleseismic  training  data  set  (parameter  numbers  from  Table  4). 


Fig.  32.  Examples  of  typical  teleseismic  records  from  the  DSI  station  of  the  ISN;  deep 
earthquake,  nuclear  test  and  shallow  earthquake. 
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For  calculating  the  P/P-coda  parameter  the  only  automatically  selected  channels  are  those 
with  SNR>3.  As  a  result  there  were  3  to  30ISN  channels  per  event  which  were  used  to 
calculate  the  statistic  average.  The  above  mentioned  scalar  discriminators  and  their 
combinations  are  presented  in  Table  4. 

Results  of  the  teleseismic  events  classification  are  presented  in  Appendix  Tables  A3  and 
A4  (together  with  depth  specification  and  empirical  thresholds),  Table  4  and  Fig.  3 1  C,D.  From 
the  3-voice  majority  voting  analysis  it  follows  that: 

•  Nuclear  tests  have  on  average  larger  values  of  P/P-coda  ratio  than  earthquakes  and  lower 
values  of  Semblance  (0.6-2  Hz)  and  Energy  Spectral  ratio  (0.6-1  Hz)/(l-3  Hz),  but  no  single 
parameter  separates  the  data  sets  completely. 

•  Each  threshold  corresponds  to  a  pair:  number  of  wrongly  classified  nuclear  tests  and  number 
of  wrongly  classified  earthquakes.  For  empirically  selected  thresholds  these  numbers  are: 
{0,9}  for  P/P-coda  ratio,  {5,5.5}  for  Semblance  1,  {1,6}  for  Spectral  ratio  1,  {2,1}  for 
Semblance  2  and  {6,4}  for  Spectral  ratio  2. 

•  When  only  shallow  earthquakes  are  considered  the  relative  number  of  errors  is  reduced: 
{0,2}  for  P/P-coda  ratio,  {5,2.5}  for  Semblance  1,  {1,2}  for  Spectral  ratio  1,  {2,0}  for 
Semblance  2,  {6,3}  for  Spectral  ratio  2. 

•  Voting  procedure  for  the  case  “with  prefiltering”,  applied  to  the  whole  data  set,  results  in 
five  errors:  earthquakes  QTIO  (deep  Afghanistan),  QT30  (deep  Tibet),  QT25  and  QT38 
(shallow,  China)  and  QT36  (shallow,  Kirgizia)  with  no  misclassifications  for  nuclear  tests. 
If  only  27  shallow  earthquakes  are  considered,  we  observe  identification  errors  for  three 
earthquakes. 

•  Voting  procedure  for  the  case  “without  prefiltering”,  applied  to  the  whole  data  set,  results  in 
two  errors:  nuclear  test  NT18  (Novaya  Zemlya)  and  earthquake  QT30  (deep  Tibet). 

•  If  only  shallow  earthquakes  are  used  for  classification,  then  the  combination  of  all  five 
parameters  mentioned  is  error  free  as  a  result  of  majority  voting  (see  Table  4). 

The  LDF  and  ANN  classifiers  applied  to  the  three  selected  scalar  discriminators  and  their 
combinations  were  tested  on  the  reduced  data  set  comprising  29  nuclear  tests  and  27  shallow 
earthquakes.  The  results  of  the  leave-one-out  cross-validation  test,  presented  in  Table  4  and 
Fig.  3 1  C,D  are  as  follows: 

•  Each  of  the  scalar  discriminators  provides  some  misclassifications  for  {LDF,  ANN}  on  the 
selected  data  base:  P/P-cof/a  -  {3,5}  mistakes.  Spectral  mt/o/ -  {6,6}  mistakes.  Spectral 
ratio!  -  {6,10}  mistakes,  Semblancel  -  {10,8}  mistakes.  Semblance!  -  {8,5}  mistakes. 
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•  Combinations  of  the  scalar  discriminants  reduce  classification  mistakes  with  the  best  variant 
giving  one  classification  mistake:  the  vector  (T1,T2,T3,T5)+LDF.  The  error  is  in  identifying 
the  nuclear  test  NTl  8  (Novaya  Zemlya). 

Table  4.  Discrimination  results  for  the  teleseismic  dataset  (27  shallow  earthquakes  and  29 


nuclear  tests)  by  multivariate  procedures. 


Parameter 

Name,  Spectral  Bands 

Threshold 

Number  of  Errors 

LDF 

•ANN:-:.. 

T1 

SEMBLANCE  1 

Spectral  semblance  (0.6-2  Hz) 
with  prefiltering  (0.5-5  Hz) 

7.5 

10 

8 

T2 

SEMBLANCE  2 

Spectral  semblance  (0.6-2  Hz) 
without  prefiltering 

0.8 

2 

8 

5 

T3 

LOGSPECTRAL  RATIO  1: 

Logarithm  of  spectral  ratio 

(0.6-lHz)/l-3Hz) 

with  prefiltering  (0.5-5  Hz) 

-0.1 1 

3 

6 

6 

T4 

LOGSPECTRAL  RATIO  2: 

Logarithm  of  spectral  ratio 
(0.6-lHz)/l-3Hz) 
without  prefiltering 

0.47 

9 

6 

10 

T5 

P/P-CODA  RATIO 

with  prefiltering  (0.5-5  Hz) 

2.3 

2 

3 

5 

T6 

Parameter  vector  (T2,T3,T5) 

- 

3 

3 

6 

T7 

Parameter  vector  (T2,T4,T5) 

HIH 

3 

4 

T8 

Parameter  vector  (T1,T3,T5) 

- 

2 

4 

T9 

Parameter  vector 
(T1,T2,T3,T5) 

- 

3 

i 

3 

TIO 

Parameter  vector 
(T1,T2,T3,T4,T5) 

- 

0 

I 

4 

As  a  result  of  the  analysis  facilitating  the  multistation  Integrative  Approach,  the  Fisher 
Linear  Discriminator  Function  and  the  Artificial  Neural  Network  procedure,  we  selected  the 
following  features  most  informative  for  classification  of  man-made  and  natural  seismic  events 
and  drew  the  following  conclusions: 
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Local  events 


•  The  best  classification  results  with  an  error  probability  equal  to  0.006  were  achieved  by 
ANN  at  the  5  band  Spectrogram  SI . 

•  The  scalar  statistics  of  Spectral  ratio  (1-3  Hz)/(6-8  Hz)  and  the  Spectrogram  SI  (7-9  Hz), 
which  greatly  resemble  each  other,  are  the  closest  competitors.  Their  success  is  explained 
by  the  physical  phenomenon  of  the  excitation  of  the  long  period  surface  waves  by 
explosions  in  the  region  with  thick  subsurface  sediments. 

•  The  more  complicated  and  interesting  spectral  semblance  and  velogram  parameter  C 
discriminators  lag  slightly  behind  with  an  error  probability  equal  to  -0.04. 

•  The  only  event  of  145  which  is  misclassified  in  the  best  case  mentioned  above,  is  the 
Dead  Sea  earthquake  which  is  estimated  as  being  located  at  a  very  shallow  depth. 
Misclassification  occurred  due  to  the  excitation  of  low  frequency  (about  0.6  Hz)  surface 
wave  energy  observed  at  some  of  the  nearby  stations. 

•  Although  the  data  belong  to  three  different  geological  areas,  it  was  handled  without 
splitting  and  with  good  classification  results,  thus  showing  robustness  and 
transportability  of  the  classification  procedures  designed. 

Teleseismic  eventx 

•  Events  of  magnitude  5  and  higher  provide  sufficient  signal-to-noise  ratio  at  the  ISN 
stations  to  perform  reliable  classification  of  earthquakes  and  nuclear  tests  from  Eurasia. 

•  P/P-coda  ratio  classifier,  as  well  as  the  energy  spectral  ratio  seem  to  be  favorable  at 
NORESS,  appear  to  be  useful  at  ISN,  when  the  statistics,  evaluated  for  the  individual 
channels,  are  averaged  over  subnetworks  of  stations. 

•  The  new  spectral  semblance  statistic,  measuring  coherency  of  spectra  at  different 
stations,  demonstrates  high  performance  in  the  classification  of  earthquakes  and  nuclear 
tests  in  the  low  frequency  range  (0.6-2  Hz).  The  effect  may  be  explained  by  frequency 
dependent  subsurface  scattering. 

•  BP  filtering  of  the  seismograms  (0.5-5  Hz)  enhances  performance  of  the  spectral  ratio 
statistic  and  reduces  performance  of  the  semblance  statistic. 

•  Joint  application  of  the  five  selected  parameters  in  the  majority  voting  procedure 
provides  full  separation  of  the  data  set,  comprising  nuclear  tests  and  shallow  earthquakes. 
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6.  APPLICATION  OF  ARRAY  METHODS  TO  TELESEISMIC  ISN 
RECORDS. 

6.1  Characterization  of  the  Array  Technique. 

Array  technique,  namely  beamforming  and  F-K  analysis  applied  to  the  closely  spaced 
and  synchronized  stations  has  the  potential  to  extract  signals  with  greater  detectability  than  the 
network  technique.  This  has  been  proved  by  the  experience  obtained  from  the  European 
seismic  arrays  such  as  NORESS,  ARCESS  in  Norway  (Mykkeltveit  et  al.,  1990),  FINESS  in 
Finland  (Uski,  1990),  GERESS  in  Germany  (lost,  1996)  and  others. 

Two  main  conditions  are  required  for  realization  of  this  advantage:  coherency  of  signals 
at  different  stations  and  stationary  noise.  Statistically,  in  the  case  when  signal  is  coherent  and 
noise  is  stationary  multivariate  random  process  with  independent  equal  variance  components, 
the  SNR  is  enhanced  at  the  output  of  a  beamformer  as  a  factor  of  N,  where  N  is  the  number  of 
stations.  Experimentally  the  standard  beamforming  has  proved  to  be  efficient  at  the  NORESS- 
type  microarrays  (25  short  period  stations,  ~3  km  radius)  for  enhancing  weak  event  signals. 
Utilization  of  the  noise  covariance  matrix  function  leads  to  an  increase  in  microarray 
efficiency.  As  shown  by  Kushnir  etal.  (1990),  application  of  the  adaptive  beamforming 
technique  (ABT)  to  the  microarray  data  may  enhance  SNR  dramatically  (much  more  than  a 
factor  of  N),  provided  the  noise  is  coherent. 

The  large  aperture  arrays,  such  as  NORSAR  (100*100  km  area)  are  probably  less 
effective  in  signal  detection,  but,  for  relatively  large  events,  provide  far  greater  precision  in 
location  than  microarrays.  Similar  results  are  expected  from  the  application  of  the  array 
procedures  to  ISN  subnetworks.  This  is  explained  by  the  fact  that  the  signal  is  no  more 
coherent  at  the  stations  separated  by  tens  of  kilometers  even  for  teleseismic  events  and  noise  is 
often  not  stationary  and  not  coherent  at  the  time  intervals  of  interest. 

An  important  feature  of  the  regional  network  is  that  the  level  of  noise  at  different  stations 
may  vary  considerably,  so  that  it  may  be  reasonable  to  expect  higher  detection  capabilities 
from  the  best  station  rather  than  by  applying  the  array  technique.  We  applied  array  technique  to 
the  ISN  channels  and  carried  out  the  following  experiments: 

•  Investigation  of  pure  noise  suppression  and  weak  signals  extraction  using  beamforming 
and  adaptive  beamforming  techniques  for  the  ISN  subnetworks. 

•  Apparent  velocity  and  azimuth  estimation  for  the  relatively  small  teleseismic  events 
using  advanced  statistical  methods  for  the  P-phase  onset  estimation  and  alignment  of 
records. 
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6.2  Regional  Network  Noise  Investigation 

The  above  question  can  be  statistically  formulated  in  terms  of  the  following: 

Theorem:  Let  us  assume  that  noise  components  are  independent  stationary  Gaussian  time 
series.  Let  W,  denote  spectral  density  N-number  of  stations)  at  a  fixed  frequency,  and 

the  signal  is  a  scalar  Gaussian  stationary  process  (true  approximation  in  a  small  time  window) 
similar  for  all  channels  S(t)=Sj  (t),  'with  spectral  density  F  at  the  same  frequency  as  above.  Let 
us  evaluate  SNR  at  the  output  of  beamformer  as  spectral  ratio.  Then: 

SNR  =  N'F/E(Wj)<NF/W„i„  (6.1) 

and  SNR  at  the  output  of  adaptive  beamformer  will  be 


SNR  =  FE(lAVj)>F/W„|„.  (6.2) 

Corollary:  If  NW^jj^  is  less  than  average  Wj,  then  it  is  better  to  select  a  channel  'with  a 
noise  level  lower  than  the  beam,  but  an  adaptive  beamformer  will  always  be  better.  From 
Eqs.  (6.1)  and  (6.2)  it  is  may  also  be  seen  that  if  the  noise  level  is  the  same  and  uncorrelated  at 
all  the  stations,  the  beam  SNR  =  F*NAV  ,  is  equal  to  the  SNR  at  the  output  of  the  adaptive 
beamformer. 

Experimental  results.  We  selected  a  2.6  minute  interval  of  pure  noise  from  14  stations 
spread  over  the  ISN  as  shovm  on  Figs.  33a  and  33b,  where  traces  are  scaled  together.  A 
stationary  noise  character  in  most  of  the  traces  and  a  lack  of  correlation  between  stations  are 
observed  and  from  Fig.  33b  we  see  that  the  noise  level  from  station  to  station  differs 
significantly.  The  standard  beamforming  technique  (SET)  of  traces  (simple  stack)  and  adaptive 
beamforming  technique  (ABT)  were  applied  to  the  different  ISN  subnetworks. 

The  results  of  application  of  the  SET  and  AET  (10  coefficients,  100  sec  adaptation 
interval)  to  all  the  14  stations  are  shown  on  Fig.  34.  Station  HRI  has  the  lowest  noise  level 
(three  times  less  than  the  second  to  it,  PRNI  station).  From  Fig.  34a  it  is  obvious  that  the  SET 
output  noise  level  is  worse  than  that  for  the  best  HRI  channel,  while  the  AET  has 
approximately  the  same  noise  level  as  HRI.  When  HRI  is  excluded  from  observations,  PRNI 
becomes  the  channel  with  minimum  noise  level,  but  there  are  nine  channels  (3,  6-14)  with 
approximately  the  same  noise  level  (Fig.  33b). 


/y\/y|fl/y^.vV^^-v^■JyW^^^ 


CBSSiKl 


VWV^/^/iA/^^VlWW 


\AV»\(VM|ANW|/Wl/^pf*»/W^ 


^■V|^'■n^vM/<''^/'/'/V^'^''^/”ii/v/'''VvA^^VA^'Wv\/^y''^y^ 


^/l%viV'/'^■*#•J^/'^/^rf\^Vl/'AwY\^^/^l\'y|^V/V'*/’Mi^'v■^'\v^/V^^^ 


■'/Iaa^'NV\/^W1|Vva^'"A/v^/‘/vY'A^'''’'*^\'''‘-'''*'^^^ 


Fig.  33.  Seismic  noise  recordings  at  the  14  ISN  stations  before  (a),  and  after  (b) 
normalization. 


Fig.  34.  Application  of  the  standard  beamforming  (SBT)  and  the  adaptive  beamforming 
(ABT)  to  ISN  noise  for:  (a)  14  channels  including  the  most  “quiet”  HRI  station;  (b) 
13  channels  without  HRI;  (c)  9  channels  with  approximately  equal  noise  level. 


Fig.  34b  illustrates  application  of  the  SBT  and  ABT  to  the  remaining  ISN  chaimels  (13). 
We  see  that  the  beam  trace  here  is  still  worse  than  the  channel  with  minimum  noise  level  and 
the  ABT  is  the  best  trace.  While  applied  to  the  nine  channels  (3,  6-14)  with  approximately 
equal  noise  level,  the  beam  and  the  ABT  have  approximately  equal  noise  levels  at  the  output, 
but  less  than  that  at  PRNI  station,  having  minimum  noise  level  in  the  selected  subnetwork  (see 
Fig.  34c) 

The  experimental  results  agree  with  the  theoretical  expectations  that  in  the  case  where 
noise  levels  vary  significantly  from  station  to  station,  selection  of  the  station  with  minimum 
noise  level  should  be  preferred  to  the  evaluation  of  beam,  but  an  adaptive  beamformer  will 
provide  results  which  are  not  worth  than  simple  selection  of  the  best  channel. 

We  tested  also  detectability  of  the  array  methods  using  ISN  records  of  real  signals  of 
relatively  weak  and  remote  events.  On  the  Fig.  35  there  is  an  example  of  the  seven  ISN  stations 
records  of  the  small  (mb=5.4)  teleseismic  event  QT18  from  Afghanistan.  The  simple  waveform 
and  high  coherency  of  the  signal  may  be  observed  here.  The  SNR  for  the  individual  ISN 
channels  is  low  and  was  enhanced  by  the  application  of  the  array  technique  (see  Fig.  35b), 
yielding  a  higher  SNR  than  for  that  of  the  best  channel  (station  RMN).  The  ABT  (25  sec  of 
adaptation,  5  coefficients)  further  improves  the  SNR.  The  traces  were  preliminary  aligned 
using  the  method  described  below.  For  another  event  QT25  (mb=5.3)  from  China,  the  results  of 
application  of  the  beam  and  ABT  are  similar  to  the  previous  example  (Fig.36). 

6.3  Apparent  Velocity  and  Azimuth  Evaluation. 

The  beam  detector  is  based  on  scanning  the  SNR  for  arriving  seismic  phases  over  a  pre¬ 
specified  grid  of  slowness  S=(Sx,Sy).  The  beam  largest  SNR  provides  estimate 
corresponding  seismic  phase  apparent  velocity  v=l/|S„axi  aJtid  azimuth  Az(Sn,ax)  to  the  source. 
Each  grid  point  Sx=|Sl*cos(Az),  Sy=|S|*sin(Az)  corresponds  to  the  specific  delay  of  signal 
arrivals  to  the  array  stations  due  to  their  different  locations  in  space.  Thus  the  estimator  can  be 
found  through  a  grid  search  to  which  best  frequency  band  may  be  added,  through  prefiltering 
by  a  series  of  band-path  filters  or  Fourier  transform.  The  large  scale  scanning  renders  this 
search  rather  exhausting  and  complicated  for  providing  accurate  solutions,  which  is  why  we 
preferred  to  use  somewhat  different  logic. 

For  the  alignment  of  seismograms  we  used  an  automatic  procedure  comprising  high 
sensitivity,  1 -component  optimal  detector  (OD)  (Kushnir  etal,  1983),  accurate  1 -component 
maxirmim  likelihood  on-set  estimator  (MLOE)  (Pisarenko  et  al.,  1987)  and  advanced  phase 
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Fig.  35.  Aligned  traces  of  the  small  (mb=5.4)  teleseismic  event  QT18  from  Afghanistan 
at  7 ISN  stations  (a);  appHcation  of  standard  beam  (SBT)  and  adaptive  beamforming 
techniques  (ABT)  to  the  ISN  subnetwork  and  comparison  with  the  “best”  RMN 
station  (b). 
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correlation  technique  (PCT)  (Gelchinsky  etal.,  1985),  thus  estimating  accurately  the  relative 
shift  between  P-wave  arrivals  to  different  ISN  stations.  The  OD  provides  a  preliminary 
estimate  of  the  P-wave  shift  dT  between  the  stations  by  evaluating  the  position  of  its  output 
maximum  (see  Fig.  37a).  The  estimate  is  corrected  by  applying  the  MLOE  on  the  time  interval 
between  the  trace  starting  point  and  the  evaluated  position  of  the  OD  output  maximum.  Finally^ 
a  small  uncertainty  time  interval  is  used  to  evaluate  the  maximum  Phase  Correlation  Function 
for  different  pairs  of  records.  The  final  vector  shift  dT  is  estimated  as  the  sum  of  shifts  at  all 
three  steps.  The  accuracy  of  the  alignment  is  checked  by  stacking  the  traces  (see  Fig.  37b).  The 
Apparent  velocity  and  Azimuth  are  finally  estimated  as  the  least  square  solution  fi-om  the 
system  of  linear  equations 


SR  =  dT  (6.3) 

where  S  is  the  unknown  slowness  and  R  is  the  matrix  of  the  stations  relative  local 
coordinates  if=(Xj-Xo,yj-yo).  The  results  of  the  slowness  estimation  for  two  weak  teleseismic 
events  (the  Semipalatinsk  nuclear  test  NT17,  mb=5.3  and  the  Chinese  earthquake  QT24, 
mb=5.2)  are  presented  in  Table  5  together  with  theoretical  evaluations  firom  the  standard  Earth 
model.  Records  of  the  seven  stations  were  used  in  the  first  case  and  8  stations  in  the  second. 
From  the  table  it  is  seen  that  in  spite  of  the  small  SNR  for  many  of  the  channels  the  estimations 
of  Slowness  and  Azimuth  are  very  close  to  the  theoretical  ones. 


Table  5.  Slowness  and  azimuth  estimated  from  the  P-phase  alignment. 


EVENT 

EARTH  ISdODEL 

REF.  STATION 

LEAST  SQUARE  ■  : 
SOLUTION 

No 

:  Lat. 

Long. 

Dist. 

Slow. 

sec/kiri 

Long. 

Slow. 

sec/km 

Azitn. 

RMS 

■ 

50.08 

78.98 

48.58 

ATZ 

32.82 

35.28 

0.079 

49.07 

0.14 

2 

44.12 

83.86 

34.63 

0.071 

56.46 

0.43 
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Fig.  37.  Apparent  velocity  and  azimuth  evaluation  for  the  small  Kazakh  nuclear  test 
NT  17:  (a)  application  of  the  Optimal  Detector  (OD)  and  Maximum  Likelihood 
Onset  Estimator  (MLOE)  to  ZNT  station  of  the  ISN;  the  traces  alignment  for  8 

ISN  stations  by  application  of  phase  correlation  algorithm  (the  bottom  trace  is  the 
stack).  65 


Conclusions:  Theoretical  analysis  and  experimental  results  of  array  technique  application  to 
subnetworks  of  Israel  Seismic  Network  show  us  that: 

•  Standard  beamforming  is  hardly  useful  for  enhancing  array  detection  since  it  may  be  more 
reasonable  to  accept  detections  from  those  stations  with  the  best  SNR. 

•  Adaptive  beamforming  techniques  may  be  useful  for  teleseismic  detection  because  it 
provides  SNR  at  the  output  no  worse  than  the  subnetwork  best  channel;  however, 
computational  cost  should  be  additionally  analyzed. 

•  We  tested  new  techniques  for  estimating  the  slowness  and  azimuth  of  teleseismic  events 
using  subnetworks  of  the  ISN  and  algorithms  for  detection  and  alignment  of  signal  traces  as 
well  as  phase  correlation  technique.  For  the  two  small  events  with  SNR  ~  3  and  with  known 
source  locations,  comparison  with  theoretical  results  showed  deviation  of  slowness 
estimation  <0.003  sec/km  (0.3 1  km/sec  of  apparent  velocity)  and  deviation  of  azimuth  <0.5 
deg,  corresponding  to  a  location  error  of  about  30-40  km  at  a  distance  of  3500-4000  km. 

7.  CONCLUSIONS  AND  RECOMMENDATIONS. 

The  potential  of  the  Israel  Seismic  Network  to  contribute  to  the  CTBT  in  detection, 
location  and  identification  of  clandestine  nuclear  tests  was  investigated  using  seismic  signals  at 
local  and  teleseismic  distance  ranges.  The  system  was  tested  for  different  kinds  of  local  events: 
ripple-fired  and  instantaneous  quarry  blasts,  earthquakes  and  underwater  explosions,  with 
considerable  success.  The  magnitude  range  for  regional  events  was  Ml~  (1 .0-3.1)  and  for  the 
Eurasian  nuclear  tests  and  distant  earthquakes  m^  ~  (5. 3-6.5). 

The  success  of  classification  is  attributed  mainly  to  the  fact  that  through  analysis  of 
dynamic  and  kinematic  features  of  the  mentioned  events,  we  managed  to  predict  theoretically 
and/or  observe  characteristic  physical  phenomena  related  to  the  concrete  class  of  sources  and  to 
design  some  statistics  based  on  seismograms  or  spectra  {velogram  parameter  C,  spectral 
semblance,  energy  spectral  ratio  and  spectrogram  determined  for  specific  frequency  bands, 
P/Pcoda)  that  most  accurately  measure  these  features. 

Analysis  of  the  physical  phenomena  related  to  explosive  sources  contributed  to  correct 
interpretation  of  the  estimated  statistics  values  and  results  of  discrimination  between  explosions 
and  earthquakes.  Preliminary  experimental  relationships  (charge  weight  vs  magnitude,  and 
peak  particle  velocity  vs  scaled  distance)  for  local  instantaneous  quarry  blasts  and  underwater 
Dead  Sea  explosions  (demonstrating  outstanding  high  coupling)  provide  a  basis  for  the  design 
and  execution  of  effective  calibration  explosions  currently  lacking  in  the  region. 
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For  the  local  events  the  characteristic  feature  seems  to  be  excitation  of  Rg,  Lg  surface 
waves  by  the  near  surface  explosive  sources,  which  together  with  ripple  firing  in  quarry  blasts 
and  bubbling  effect  in  the  underwater  explosions,  imply  enhancement  of  low  frequencies  in 
signal  spectra,  spectral  modulation  and  coherency  of  spectral  shapes  at  different  stations. 
Additional  coherency  of  the  local  explosions  spectra  may  be  attributed  to  the  spherical 
symmetry  of  the  source. 

Teleseismic  events  are  characterized  by  the  impulsive  P-onset  and  shifting  of  spectral 
energy  to  higher  frequencies  for  nuclear  detonations.  Application  of  semblance  revealed 
impredictable  larger  coherency  of  earthquake  spectral  shapes  compared  to  those  of  nuclear 
tests,  possibly  due  to  scattering  effects. 

Due  to  the  dense  distribution  of  the  ISN  stations  over  Israel,  we  had  an  opportunity  to 
accumulate  information  contained  in  the  ISN  records  and  perform  reliable  event  discrimination. 
Joint  analysis  of  the  above  mentioned  discriminants  and  averaging  of  feature  vectors  over  the 
ISN  sunbetworks  provides  additional  discrimination  effects  and  yields  an  almost  correct 
separation  of  both  local  and  teleseismic  data  bases.  Both  data  bases  were  large  enough  to  test 
performance  of  different  automatic  classifiers  and  select  the  best  of  them  according  to 
minimum  error  probability  estimator  criterion:  ANN  plus  averaged  spectrogram  in  the  first 
case  and  LDF  plus  specific  combination  of  the  teleseismic  discriminants  in  the  second  case. 

The  relatively  dense  character  of  the  ISN  subnetworks  is  attractive  for  the  application  of 
array  techniques  to  enhance  the  network  detectability  and  accuracy  of  location  in  far  regional 
and  teleseismic  ranges.  Possible  enhancement  of  the  ISN  detectability  may  be  achieved  by 
adaptation  to  the  current  multichannel  noise  characteristics  and  studying  the  signal  propagation 
characteristics  across  the  ISN  subnetworks  to  provide  minimum  detour  of  signal  dynamic  and 
kinematic  features.  As  a  result  of  the  investigation  of  the  ISN  pure  noise  and  weak 
seismograms  of  distant  events,  we  proved  the  advantages  of  the  adaptive  beamforming 
technique  over  standard  beamforming  for  the  detection  of  teleseismic  signals.  Attempts  to 
estimate  azimuth  and  apparent  velocity  of  weak  teleseismic  events  were  quite  successful  in  the 
cases  tested.  The  analysis  showed  that  for  the  purposes  of  slowness  vector  estimation,  a 
relatively  limited  station  aperture  should  be  used  to  avoid  influence  of  regional  crustal  and 
mantle  inhomogeneities. 

For  confirmation  of  the  results  and  implementation  of  the  automated  detectors  and 
classifiers  in  routine  data  processing,  we  recommend  that  they  be  tested  in  large  scale 
experiments  comprising  hundreds  of  events  from  different  regions. 
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Table  Al.  Underground  nuclear  tests  used  in  the  study, 


Location 

Date  hrmnsec 

Kazakh 

850615  005700.7 

Long, 

deg 


8.88 


*NT3  Ka2akh 


*NT4  Kazakh 


*NT5  Kazakh 


*NT6  China 


*NT8  Kazakh 


*NT9  Kazakh 


NTl  1  Kazakh 


NT12  Kazakh 


NTl  3  Kazakh 


NTl  5  W.Siberia 


NT 


NT 


NTl  8  Nov.Zemlya 


NTl  9  Kazakh 


NT21  Kazakh 


850720  005314.5 


870312  015717.3 


870417  010304.8 


870506  040205.5 


870605  045958.3 


870606  023706.9 


870802  005806.8 


871115  033106.7 


871213  032104.8 


871227  030504.7 


880213  030505.9 


880403  013305.8 


880504  005706.8 


880822  161958.2 


880914  035957.4 


881112  033003.7 


881204  051953.0 


881217  041806.9 


890122  035706.5 


890212  041506.7 


901024  145758.1 


49.951 

78.829 

49.932 

78.785 

49.886 

78.691 

49.777 

78.089 

41.584 

88.737 

49.865 

78.143 

49.880 

78.917 

49.871 

78.791 

NT23 

China 

920521  045957.5 

NT24 

China 

931005  015956.5 

NT25 

China 

940610  062557.8 

NT26 

China 

941007  032558.1 

NT27 

China 

950515  040558.0 

NT28 

China 

950817  010004.8 

*NT29 

China 

960608  025558.0 

49.864 


49.954 


49.928 


66.316 


49.833 


50.078 


73.387 


49.886 


49.895 


73.361 


41.604 


41.647 


41.527 


41.662 


41.665 


41.600 


78.758 


78.910 


78.769 


78.548 


78.808 


78.988 


54.998 


78.926 


78.758 


54.707 


88.813 


88.681 


88.753 


88.821 


Table  A2.  Selected  Eurasian  earthquakes. 


No. 

Location 

Date  hrmnsec 

Lat, 

deg 

Long, 

deg 

Depth, 

km 

Dist, 

deg 

Azim, 

deg 

mb 

*QT1 

Uzbek 

840319  202838.2 

40.320 

63.350 

15 

24 

61 

6.5 

841026  202221.8 

39.155 

71.328 

33 

30 

65 

6.0 

QT3 

China 

850823  124156.1 

39.431 

75.224 

7 

33 

65 

*QT4 

China 

850911  204549.5 

39.356 

75.407 

15 

33 

65 

5.8 

QT5 

Tajik 

851013  155951.2 

40.301 

69.823 

16 

29 

63 

5.8 

**QT6 

Uzbek 

860325  234934.4 

40.336 

63.658 

33 

24 

61 

5.2 

*QT7 

Afghanistan 

860426  141507.6 

36.495 

71.114 

187 

29 

71 

5.6 

*QT8 

Afghanistan 

860821  013417.3 

36.472 

71.083 

235 

29 

71 

*QT9 

Afghanistan 

860915  214229.2 

36.714 

71.092 

89 

29 

70 

5.8 

QTIO 

Afghanistan 

860917  120809.4 

37.290 

71.730 

120 

30 

69 

5.5 

*QT11 

China 

870225  195636.4 

38.029 

91.144 

33 

45 

66 

m 

*QT12 

China 

870430  051737.1 

39.770 

74.590 

9 

32 

*QT13 

China 

870918  215836.6 

47.276 

89.674 

33 

43 

53 

5.3 

*QT14 

Afghanistan 

871003  110005.2 

36.454 

71.437 

95 

30 

71 

5.9 

*QTI5 

Tajik 

880326  225842.8 

38.309 

73.234 

121 

31 

67 

5.7 

QT16 

Tajik 

880720  062051.4 

37.028 

72.914 

41 

31 

69 

5.5 

QT17 

Afghanistan 

880806  090321.9 

36.461 

71.043 

195 

29 

71 

6.1 

QT18 

Afghanistan 

880909  211236.1 

36.483 

71.381 

188 

30 

71 

QT19 

China 

880923  044640.6 

39.570 

74.506 

33 

32 

65 

QT20 

Afghanistan 

880926  071700.2 

36.294 

71.374 

107 

30 

71 

5.6 

QT21 

Afghanistan 

900205  051645.1 

37.047 

71.250 

110 

ll^ll 

6.1 

QT22 

China 

900417  015933.4 

39.436 

74.900 

33 

32 

65 

6.0 

*QT23 

Kazakh-China 

900803  091539.6 

47.949 

84.958 

19 

Ell 

51 

6.1 

QT24 

China 

901024  233815.1 

44.117 

83.856 

20 

39 

57 

5.2 

QT25 

China 

901024  234657.6 

44.119 

83.876 

22 

El 

5.3 

♦QT26 

China 

910819  060553.0 

46.954 

85.333 

44 

40 

53 

5.7 

**QT27 

China 

930202  160510.7 

42.212 

86.126 

10 

OI 

5.7 

QT28 

Afghanistan 

930918  050227.2 

36.369 

71.600 

117 

30 

71 

6.1 

QT29 

China 

931002  084232.8 

38.141 

88.638 

16 

43 

66 

6.3 

QT30 

Tibet 

860110  034629.9 

28.648 

86.527 

55 

44 

80 

Bl 

QT31 

Tibet 

860620  171246.9 

31.240 

86.847 

33 

43 

76 

m\ 

QT32 

China 

860826  094300.3 

n 

8 

m 

6.2 

China 

870810  121217.8 

38.118 

106.357 

10 

57 

62 

Bl 

N.China 

880103  213225.5 

38.111 

106.336 

14 

57 

62 

m\ 

QT35 

China 

881105  021430.3 

BBI 

8 

in 

5.9 

QT36 

Kirgiz 

881221  082103.7 

41.220 

72.298 

33 

31 

61 

Bl 

74 


QT38  China 


QT39  China 


31.545 

102.464 

33 

56 

71 

6.0 

36.206 

106.346 

10 

5.0 

37.819 

91.971 

12 

46 

66 

6.1 

32.432 

r  92.743 

13 

73 

5.6 

47.903 

84.961 

33 

40 

51 

5.0 

-  only  ISN  data  are  presented,  **  -  only  NORESS  data  are  presented. 


Table  A3.  Majority  voting  analysis  for  the  ISN  teleseismic  recordings 
with  (0.5  -  5  Hz)  prefiltering. 


Nuclear  tests 


P/Pcoda 


No 

Location 

NTI 

Kazach 

NT  2 

Kazakh 

NT  3 

Kazakh 

NT  4 

Kazakh 

NT  5 

Kazakh 

NT  6 

China 

NT  7 

Kazakh 

NTS 

Kazakh 

NT  9 

Kazakh 

NT  10 

Kazakh 

NTH 

Kazakh 

NT  12 

Kazakh 

Semblance  1 


Spec.  Ratio  1 


thr.=0.9 


Number  of  errors 


Voting 


3:0 


3:0 


3:0 


3:0 


3:0 


3:0 


3:0 


2:1 


3 


3:0 


3:0 


3:0 


3:0 


3:0 


3:0 


2:1 


2:1 


3:0 


3:0 


3 


3:0 


2:1 


3 


3 


3 


3 


3 


3:0 


3 
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Table  A3.  Continuation. 


Eartquakes 

P/Pcoda 

No. 

Location 

thr.=2.3 

QTi 

Uzbek 

15 

1.028 

+ 

QT2 

Tajik 

33 

1.915 

+ 

QT3 

China 

7 

0.802 

+ 

QT4 

China 

15 

1.255 

+ 

QT5 

Tajik 

16 

1.331 

+ 

QT6 

Uzbek 

33 

- 

- 

QT7 

Afghan 

187 

2.42 

- 

QT8 

3.74 

- 

QT9 

Afghan 

89 

1.12 

+ 

QTIO 

Afghan 

120 

3.85 

- 

QTll 

China 

33 

2.160 

+ 

QT12 

China 

9 

1.384 

+ 

QT13 

China 

33 

1.307 

+ 

QT14 

Afghan 

95 

1.89 

IBIiH 

QT15 

Tajik 

121 

4.16 

- 

41 

1.352 

+ 

QT17 

Afghan 

195 

1.76 

+ 

QT18 

Afghan 

188 

3.76 

- 

QT19 

China 

33 

1.733 

+ 

QT20 

Afghan 

107 

2.8 

- 

QT21 

Afghan 

no 

1.93 

4* 

QT22 

China 

33 

1.304 

+ 

QT  23 

Kaz 

19 

1.241 

+ 

QT24 

China 

20 

1.330 

+ 

QT  25 

China 

22 

1.188 

+ 

QT26 

China 

44 

1.763 

+ 

- 

- 

QT28 

Afghan 

117 

1.66 

+ 

QT29 

China 

16 

1.288 

+ 

QT30 

Tibet 

- 

QT3I 

Tibet 

33 

2.232 

+ 

QT32 

China 

8 

1.973 

+ 

QT  33 

China 

10 

0.543 

+ 

QT34 

N.China 

14 

1.506 

+ 

QT35 

China 

8 

1.375 

+ 

QT36 

Kirgiz 

33 

4.522 

- 

QT37 

China 

33 

1.118 

+ 

QT38 

China 

10 

3.301 

- 

QT39 

China 

12 

1.351 

QT40 

Tibet 

13 

1.108 

+ 

QT41 

Kaz-China 

33 

1.238 

+ 

Nunnber  of  errors:  shallow/all 

2/7 

7i 


2.5/3 


2/4 


Table  A4.  Majority  voting  analysis  for  the  ISN  teleseismic  recordings 
without  prefiltering. 


Nuclear  tests 

P/Pcoda 

Semblance  2 

Spectral  Ratio  2 

No. 

Location 

value 

thr.=1.6 

Voting 

NTl 

Kazach 

0 

3.844 

+ 

0.537 

4- 

0.853 

4- 

3:0 

NT  2 

Kazakh 

+ 

0.434 

+ 

0.821 

4- 

3:0 

NT  3 

Kazakh 

0 

3.781 

+ 

0.638 

4- 

0.676 

4- 

3:0 

NT  4 

Kazakh 

0 

4.260 

+ 

0.760 

+ 

1.870 

2:1 

NT  5 

Kazakh 

+ 

0.782 

4- 

0.988 

4- 

3:0 

NT  6 

China 

0 

IIQH 

+ 

0.104 

4- 

0.483 

4* 

3:0 

NT  7 

Kazakh 

+ 

0.656 

4- 

1.003 

4- 

3:0 

NT  8 

Kazakh 

0 

4.274 

+ 

0.870 

- 

1.183 

4- 

2:1 

NT  9 

Kazakh 

0 

4.302 

+ 

0.615 

4- 

1.075 

4- 

3:0 

NT  10 

Kazakh 

+ 

0.588 

4* 

1.568 

4- 

3:0 

NT  11 

Kazakh 

0 

4.253 

+ 

0.678 

4- 

1.701 

- 

2:1 

NT  12 

Kazakh 

+ 

0.770 

4- 

1.648 

2:1 

NT  13 

Kazakh 

0 

4.499 

+ 

0.681 

4- 

1.331 

4- 

3:0 

NTH 

Kazakh 

0 

4.029 

+ 

0.349 

4- 

0.895 

4- 

3:0 

NT  15 

W.Sib. 

3.763 

+ 

0.670 

4* 

0.530 

+ 

3:0 

NT  16 

Kazakh 

0 

3.317 

+ 

0.794 

4-- 

0.843 

4- 

2:1 

NT  17 

Kazakh 

0 

4.936 

+ 

0.599 

+ 

+ 

3:0 

NT  18 

N.Zemlya 

0 

2.842 

+ 

0.840 

- 

1:2 

NT  19 

Kazakh 

0 

5.430 

4* 

0.752 

+ 

1.364 

4- 

3:0 

NT  20 

Kazakh 

0 

4.280 

+ 

0.623 

4- 

- 

3:0 

NT  21 

Kazakh 

+ 

4- 

1.200 

4* 

3:0 

NT  22 

N.Zemlya 

0 

3.400 

+ 

0.763 

4- 

2.560 

- 

2:1 

NT  23 

China 

+ 

0.595 

4- 

1.103 

4- 

3:0 

NT  24 

China 

4- 

0.174 

4- 

0.606 

4- 

3:0 

NT  25 

China 

0 

4.142 

4- 

0.463 

4- 

0.696 

4- 

3:0 

NT  26 

China 

4- 

0.442 

4- 

0.489 

4- 

3:0 

NT  27 

China 

0 

3.800 

4- 

0.243 

4- 

0.498 

4- 

3:0 

NT  28 

China 

0 

3.990 

4- 

0.321 

4- 

0.584 

4- 

3:0 

NT  29 

China 

0 

3.611 

+ 

0.491 

4- 

0.676 

4- 

3:0 

Number  of  errors : 

2 

6 

Table  A4.  Continuation. 
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LIVERMORE,  CA  94551 


LAWRENCE  LIVERMORE  NAT’L  LAB 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
POBOX  808,  MS  L-221 
LIVERMORE,  CA  94551 


JAMES  LEWKOWICZ 
WESTON  GEOPHYSICAL  CORP, 
325  WEST  MAIN  STREET 
NORTHBORO,  MA  01532 


LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  1663,  MS  C335 
LOS  ALAMOS,  NM  87545 


BRIAN  MITCHELL 

DEPT  OF  EARTH  &  ATMOS.  SCIENCES 
ST.  LOUIS  UNIVERSITY 
3507  LACLEDE  AVENUE 
ST.  LOUIS,  MO  63103 


JAMES  NI 

NEW  MEXICO  STATE  UNIVERSITY 
DEPARTMENT  OF  PHYSICS 
LAS  CRUCES,  NM  88003 


JOHN  ORCUTT 

INST.  OF  GEOPH.  &  PLANETARY  PHYSICS 
UNIV.OF  CALIFORNIA,  SAN  DIEGO 
LA  JOLLA,  CA  92093 


PACIFIC  NORTHWEST  NAT’L  LAB 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  999,  MS  K6-84 
RICHLAND,  WA  99352 


KEITH  PRIESTLEY 
DEPARTMENT  OF  EARTH  SCIENCES 
UNIVERSITY  OF  CAMBRIDGE 
MADINGLEY  RISE,  MADINGLEY  ROAD 
CAMBRIDGE,  CB3  0EZUK 


LAWRENCE  LIVERMORE  NAT’L  LAB 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  808,  MS  L-208 
LIVERMORE,  CA  94551 


LAWRENCE  LIVERMORE  NAT’L  LAB 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  808,  MS  L-205 
LIVERMORE,  CA  94551 


THORNE  LAY 

UNIV.  OF  CALIFORNIA,  SANTA  CRUZ 
EARTH  SCIENCES  DEPARTMENT 
EARTH  &  MARINE  SCIENCE  BUILDING 
SANTA  CRUZ,  CA  95064 


LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
POBOX  1663,  MS  F659 
LOS  ALAMOS,  NM  87545 


GARY  MCCARTOR 
SOUTHERN  METHODIST  UNIVERSITY 
DEPARTMENT  OF  PHYSICS 
DALLAS,  TX  75275-0395 


RICHARD  MORROW 

USACDA/IVI 

320  2 1ST  STREET,  N.W. 

WASHINGTON  DC  20451 


ROBERT  NORTH 

CENTER  FOR  MONITORING  RESEARCH 
BOON.  17th  STREET,  SUITE  1450 
ARLINGTON,  VA  22209 


PACIFIC  NORTHWEST  NAT’L  LAB 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  999,  MS  K6-48 
RICHLAND,  WA  99352 


PACIFIC  NORTHWEST  NAT’L  LAB 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
PO  BOX  999,  MS  K5-12 
RICHLAND,  WA  99352 


JAY  PULLI 

BBN  SYSTEMS  AND  TECHNOLOGIES,  INC. 
1300  NORTH  17TH  STREET 
ROSSLYN,  VA  22209 
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DELAINE  REITER 
AFRLA^SOE  (SENCOM) 

29  RANDOLPH  ROAD 
HANSCOMAFB,MA  01731-3010 


^  DAVID  RUSSELL 
HQ  AFTAC/TTR 
1030  SOUTH  HIGHWAY  A1 A 
PATRICK  AFB,  FL  32925-3002 


SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  93 11 

MS  11 59,  POBOX  5800 
ALBUQUERQUE,  NM  87185-1159 


THOMAS  SERENO,  JR. 

SAIC 

10260  CAMPUS  POINT  DRIVE 
SAN  DIEGO,  CA  92121 


MATTHEW  SIBOL 
ENSCO,  INC. 

445  PINEDA  CT. 
MELBOURNE,  FL  32940 


BRIAN  SULLIVAN 
BOSTON  COLLEGE 
INSITUTE  FOR  SPACE  RESEARCH 
140  COMMONWEALTH  AVENUE 
CHESTNUT  HILL,  MA  02167 


LAWRENCE  TURNBULL 

ACIS 

DCI/ACIS 

WASHINGTON  DC  20505 


TERRY  WALLACE 
UNIVERSITY  OF  ARIZONA 
DEPARTMENT  OF  GEOSCIENCES 
BUILDING  #77 
TUCSON,  AZ  85721 


RU  SHAN  WU 

^  UNIV.  OF  CALIFORNIA  SANTA  CRUZ 
EARTH  SCIENCES  DEPT. 

1156  HIGH  STREET 
SANTA  CRUZ,  CA  95064 

» 


PAUL  RICHARDS 
COLUMBIA  UNIVERSITY 
LAMONT-DOHERTY  EARTH  OBSERV. 
PALISADES,  NY  10964 


CHANDAN  SAIKIA 

WOOODWARD-CLYDE  FED.  SERVICES 
566  EL  DORADO  ST.,  SUITE  100 
PASADENA,  CA  91 101-2560 


SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  5704 

MS  0655,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0655 


AVI  SHAPIRA 
SEISMOLOGY  DIVISION 
IPRG 

P.O.B.  2286 

NOLON  58122  ISRAEL 


DAVID  SIMPSON 
IRIS 

1200  NEW  YORK  AVE„  NW 
SUITE  800 

WASHINGTON  DC  20005 


TACTEC 

BATTELLE  MEMORIAL  INSTITUTE 
505  KING  AVENUE 

COLUMBUS,  OH  43201  (FINAL  REPORT) 


GREG  VAN  DER  VINK 
IRIS 

1200  NEW  YORK  AVE.,  NW 
SUITE  800 

WASHINGTON  DC  20005 


JILL  WARREN 

LOS  ALAMOS  NATIONAL  LABORATORY 
GROUP  NIS-8 
P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545  (5  COPIES) 


JIAKANG  XIE 
COLUMBIA  UNIVERSITY 
LAMONT  DOHERTY  EARTH  OBSERV. 
ROUTE  9W 

PALISADES,  NY  10964 


MICHAEL  RITZWOLLER 
DEPARTMENT  OF  PHYSICS 
UNIVERSITY  OF  COLORADO 
CAMPUS  BOX  390 
BOULDER,  CO  80309-0309 


SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  5704 

MS  0979,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0979 


SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  5736 

MS  0655,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0655 


ROBERT  SHUMWAY 
410  MRAK  HALL 
DIVISION  OF  STATISTICS 
UNIVERSITY  OF  CALIFORNIA 
DAVIS,  CA  95616-8671 


JEFFRY  STEVENS 
MAXWELL  TECHNOLOGIES 
8888  BALBOA  AVE. 

SAN  DIEGO,  CA  92123-1506 


NAFI  TOKSOZ 

EARTH  RESOURCES  LABORATORY,  M.I.T. 
42  CARLTON  STREET,  E34-440 
CAMBRIDGE,  MA  02142 


FRANK  VERNON 

UNIV.OF  CALIFORNIA,  SAN  DIEGO 
SCRIPTS  INST.  OF  OCEANOGRAPHY 
IGPP,  0225 

9500  GILMAN  DRIVE 
LA  JOLLA,  CA  92093-0225 

DANIEL  WEILL 
NSF 

EAR-785 

4201  WILSON  BLVD.,  ROOM  785 
ARLINGTON,  VA  22230 


JAMES  E.  ZOLLWEG 
BOISE  STATE  UNIVERSITY 
GEOSCIENCES  DEPT. 

1910  UNIVERSITY  DRIVE 
BOISE,  ID  83725 
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